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Abstract 

Abstract: We study a lattice field model which qualitatively reflects 
the phenomenon of Anderson localization and delocalization for real sym- 
metric band matrices. In this statistical mechanics model, the field takes 
values in a supermanifold based on the hyperbolic plane. Correlations in 
this model may be described in terms of a random walk in a highly corre- 
lated random environment. We prove that in three or more dimensions the 
model has a 'diffusive' phase at low temperatures. Localization is expected 
at high temperatures. Our analysis uses estimates on non-uniformly ellip- 
tic Green's functions and a family of Ward identities coming from internal 
supersymmetry. 



1 Introduction 

1.1 Some history and motivation 

It has been known since the pioneering work of Wegner [20l [19] that informa- 
tion about the spectral and transport properties of random band matrices and 
random Schrodinger operators can be inferred from the correlation functions of 
statistical mechanical models of a certain kind. These models have a hyperbolic 
symmetry, typically a noncompact group such as 0{p,q) or U(p, g), and were 
originally studied in the limit oi p — q — Q replicas. 
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The connection between random Schrodinger operators and statistical me- 
chanics models was made more precise by Efetov [7], who introduced the so- 
called supersymmetry method to avoid the use of replicas. In Efetov's for- 
mulation one employs both commuting (or bosonic) and anticommuting (or 
Grassmann) integration variables, and these are related by a natural symmetry 
that makes the emerging statistical mechanics system supersymmetric (SUSY) . 
The simplest class of these models has a U(l, 1|2) symmetry. This means that 
for the bosonic variables there exists a hyperbolic symmetry U(l, 1) preserving 
an indefinite Hermitian form on C^, and the Grassmann variables are governed 
by a compact U(2) symmetry. Moreover, there exist odd symmetries mixing 
Grassmann and bosonic variables. 

The fields Z'' 9 j i— > Qj of the supersymmetric models introduced by Efetov 
are 4 by 4 supermatrices built from bosonic as well as Grassmann entries. In the 
physics literature one usually assumes the sigma model approximation, which is 
believed to capture the essential features of the energy correlations and transport 
properties of the underlying quantum system. The sigma model approximation 
constrains the matrix field Q hy = Id for all j. This constraint is similar to 
the constraints appearing in the Ising or Heisenberg models, where Sj ■ Sj ~ 1 . 
We refer the reader to [TH [HI [31 [H] for an introduction to these ideas. 

The models described above are difficult to analyse with mathematical rigor 
in more than one dimension. In this paper we study a simpler SUSY model. 
Our exposition will be essentially self-contained and the full supersymmetric 
formalism alluded to here will serve primarily as a source of motivation. 

1.2 Probabilistic representation of our model 

In this paper we analyze a lattice field model which may be thought of as a 
simplified version of one of Efetov's nonlinear sigma models. More precisely, it 
is related to the model that derives from real symmetric matrices, see Section|3l 
In this statistical mechanics model the field at site j has four degrees of freedom. 
Two of these, tj and Sj , parametrize a hyperboloid and the other two, ipj and 
tjjj , are Grassmann (i.e., anticommuting) variables. Technically speaking, the 
field takes values in a target space denoted by H^'^, which is a supermanifold 
extension of the hyperbolic plane H^; see Section [51 This model was introduced 
by one of us in [21i j5j. and localization was established in one dimension (ID) in 
the sense that the conductance was proven to decay exponentially in the system 
size [21]. The model is expected to refiect the qualitative behavior of random 
band matrices - namely localization and diffusion - in any dimension. However, 
most of our discussion will be restricted to three dimensions. 

Our supersymmetric hyperbolic nonlinear sigma model, called the H^'^ model 
for short, will be formulated on a lattice cube A C Z'' of side L . We shall see 
(in Section [2T2]) that the action of the field variables is quadratic in ip, -0, and 
s. This feature is special to the horospherical coordinate system that we use. 
It enables us to reduce the H^'^ model to the statistical mechanics of a single 
field t : A ^ M. , j 1-^ tj . Its free energy or effective action, F{t), is real, so the 
resulting statistical mechanical model has a probabilistic interpretation. 
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In order to specify F{t), first consider the finite-difference elliptic operator 
Dp^i^it) defined by the quadratic form 

[v ■ D^Jt) «]a = p Y.^,,^ - «^)' + ' EfceA ^'""^l ■ (1-1) 

This operator plays a central role in our analysis. The first sum is over nearest 
neighbor pairs in A, and [; ]a denotes the usual scalar product in i^{A). We 
see that Difi{0) is the finite-difference Laplacian. The regularization parameter 
£ > will serve to make the theory well-defined. One may interpret Dp^e{t) as 
the generator of a random walk in an environment given by the fiuctuating field 
t, with a death rate of ee*^ at site j. Note that the operator D is elliptic but 
not uniformly so, as tj G M has unbounded range. 

The free energy or effective action Fp_^{t) is now expressed by 

- InDet^/^D/j^eW + ^^.^j^i^k -s + e coshtfc) . (1.2) 
If dfj,\ denotes the product measure 

= n # (1-3) 

on R'^I , then the partition function is 

ZA(/3,e) = / e-^^-dfiA = 1 . (1.4) 

The partition function is identically equal to unity independent of /?, e even 
when P depends on the edge (ij) and e depends on the lattice point k ; see 
(j5.1|) . This is a refiection of an internal supersymmetry which will be explained 
in later sections. There exist many variants of this identity. One of them gives 
us easy control of nearest neighbor fluctuations of the field t (cf. Section [6]). 
The expectation of an observable function t ^ f{t) is defined by 

{f)K,p,e = j fe-^^-^^dfij, . (1.5) 

Let us make a few comments on these expository definitions. 

1. The action or free energy F/3,e(i) is nonlocal due to the presence of the 
term - lnDet^/^L'(t). This nonlocality arises from integrating out three 
massless free fields, one (s) of bosonic and two (-i/i, -0) of Grassmann type. 

2. Fp^i^lt) is not convex as a function of t and therefore the Brascamp-Lieb 
estimates used in earlier work on a related model [18] do not apply. The 
lack of convexity is an important feature and opens the possibility for a 
localization-delocalization transition to occur. 
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3. When e = 0, F^,o(0 invariant under shifts tj tj + c by any constant 
c G M . To see this, note that for e = we have Di3^o{t + c) — e^'^£)/3,o(*) by 
(fTT]) . The resulting additional term -|A|c from - lnDet^/^D/3,o(<) in ^L2l 
is canceled by another such term, which arises from shifting J^keA ■ This 
symmetry (which is a formal one, since the integral is ill-defined for e = 0) 
is associated with the presence of a massless mode. The importance of the 
regularization e, which was omitted from the present argument, becomes 
evident from the saddle point discussed below. 

4. The model at hand describes a disordered quantum system at zero tem- 
perature. Nevertheless, adopting the familiar language of statistical me- 
chanics and thermodynamics, we refer to the field stiffness /3 as the inverse 
'temperature'. (/3 is actually the dimensionless conductance for an Ohmic 
system of size L = 1 as measured in lattice units.) 

1.3 Main result 

The main goal of this paper is to estimate the fluctuations of the field t for 
large values of the parameter f3 and dimension d = 3. This will enable us to 
prove that the random walk in the random environment drawn from F(t) is 
transient. More precisely, we will prove the following. (Similar estimates hold 
for all dimensions d > 3.) 

Theorem 1 For d — 3, there is a (3 > \ such that if (3 > (3, the fluctuations of 
the field t are uniformly hounded in x, y, and A : 

(cOsh™(t,-tj,))A,/3,e < 2 , (1.6) 

provided that m < [3^/^ . 

This theorem implies that for any x and y, \tx — ty \ is very unlikely to be large. 
A stronger version of ()1.6p is given in (19. 3|) . We will use this result to prove 

Theorem 2 Under the hypothesis of Theorem]^ the average field is bounded: 

{coshP{t,))^jj,,<^, (1.7) 

providedp < 10 and |A|^^"/'^e > 1 witha > l/ln/3. Thus in the thermodynamic 
limit |A| —^oowe may send £ — > while maintaining the bound on {cosW t^) ■ 

To investigate the localized or extended nature of the energy eigenstates of 
a disordered quantum system with Hamiltonian H, one looks at the average 
square of the quantum Green's function, \{H — E + ie)~^(x, ?/)p. The analog of 
this Green's function in the H^'^ model is the two-point correlation function 

C^y = {e*^ So: e^'' sy) , (1.8) 
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where the expectation is given by the full functional integral defined in Sections 
12. 1[ 12.21 After integration over the fields -0, ip, and s, we have 

C,, = {e'-+*yDp.At)-'{x,y))^^^^^^ ^ y))^^^^^ , (1.9) 

where D = e~*D o e~*. Note that Cxy is positive both pointwise and as a 
quadratic form. A simple calculation shows that 

Df3,,it)^-pA + f3V{t)+ee-' , (1.10) 

where V{t) is a diagonal matrix (or 'potential') given by 

(sum over nearest neighbors) and e~* is the diagonal matrix with {c^*)jj = e~*J . 
In AppendixiBjwe establish the sum rule e X^yeA ~ 1' reflecting conservation 
of probability for the quantum dynamics generated by a Hamiltonian H . 

Note that if t were bounded, then D(t) (given by (jl.ip 'l would be uniformly 
elliptic and we could establish good diffusive bounds on the two-point function 
C (|1.9p . However, Theorems [T] and [2] only say that large field values are unlikely. 
To get optimal bounds on C we would need to prove uniform ellipticity on a 
percolating set. The set on which \tj + tji \ < M, is presumably a percolating 
set but this does not readily follow from our estimates. 

Our next theorem states a quasi-diffusive estimate on C. More precisely let 
Go = {—f3A + s)~^ be the Green's function for the discrete Laplacian (with a 
regularization term e) and Go = (— /3A -I- e/2)~^. In 3 dimensions Gs{x,y) < 
/3~^(1 + \x ~ y\)^^ (and the same is true for Go). Then we have 

Theorem 3 Let / : A ^ M fee non-negative. Then assuming the hypotheses of 
Theorems\^ and\^ we have 

^[f; Go/] < [/; Cf] = ^ a, /(z)/(j) < K[f; Go/] , (1.11) 

where f{j) = (1 + \j — a;|")^^/(j), x CI A is any fixed point, and K and K' are 
constants independent of f . 

Remark. In this paper we always use periodic boundary conditions on A C 
1? . The distance |x — y| between two points is always the distance on A with 
periodic boundary conditions. 

1.4 Saddle point 

One may try to gain a crude understanding of the behavior of the H^'^ sigma 
model via a simple saddle-point analysis. Let be the configuration of < = {tj} 
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which minimizes the effective action Fp^^{t) defined in (|1.2p . In Appendix [XI we 
prove that is unique and t^^^ = t* independent of j . For large (3 we find 

ID: ee-**~/3-\ 2D: ee-^'-e"'', (1.12) 

in one and two dimensions, respectively. Thus in ID or 2D the saddle point 
depends sensitively on the regularization parameter e. The value of t* suggests 
a strong asymmetry of the field favoring negative values of t. On the other hand, 
in 3D at low temperatures, we find t* — Q independent of e. Our estimates ()1.7p 
confirm this value by controlling fluctuations about the saddle. For (3 small, in 
3D, the saddle t* is again strongly e-sensitive, suggesting localization. 

The bias to negative values of the field t is expected to be closely related 
to localization. Note that since —A + V{t) > 0, the additional term £e~* 
makes Dp^i. strictly positive at the saddle suggesting that C^y decays roughly 
like e~™l^~^l with m? = ee~* / (3 = /3~^ and in ID and 2D respectively. 
There are important fluctuations away from this saddle but we do not expect 
them to spoil the exponential decay. For the ID chain this has been proved [21]. 

1.5 Edge reinforced random walk 

A number of mathematicians (Kozma, Heydenreich, Sznitmann) have noted 
that our random walk looks similar to a linearly edge reinforced random walk 
(ERRW) . ERRW is a history-dependent walk which prefers to visit edges it has 
visited in the past. Let n(e) denote the number of times the walk has visited the 
edge e. Then the probability that the walk at vertex v will visit a neighboring 
edge e equals {a+n{e)) / Sa{v) where S is the sum of a+n(e') over all the edges e' 
touching V. The parameter a is analogous to our (3. Coppersmith and Diaconis 
[2] proved that this history-dependent walk can be expressed as a random walk 
in a random environment; see also more recent work by Merkl and RoUes |13j 
in which recurrence of the walk is established on a 2D lattice for small (3. This 
is analogous to localization in our model. The environment of ERRW is very 
similar to the environment in H^'^. In fact, both environments have nonlocal 
actions arising from the square root of a determinant. Although the two models 
do not seem to be identical, they may have similar properties. 

1.6 Outline of the paper 

The remainder of this paper is organized as follows. In the next section we 
give a precise definition of the full H^'^ model and introduce the horospherical 
coordinate system. The effective action defined in (|1.2|) is then derived by 
integration of the field s and the Grassmann fields t/j and ■(/' • Section [3] provides 
a physical motivation for the study of this model. In Section |4] we explain the 
symmetries of the model and briefly discuss its perturbative renormalization 
group flow. The basic Ward identities we shall need are given in Section [5] and 
are derived in Appendix [C] Section [5] ends with a rough outline of our proof 
and a description of the remaining sections of this paper. 
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2 Definition of the model 

We now fill in the details of the definition the H^'^ model and derive the free 
energy Fj^^i^it) given above. 

2.1 Full supersymmetric model 

As in Section [L2l let A C Z'* be a cube of size L. For each lattice site j g A we 
introduce a supervector Uj g R'^l^, 

= izj,Xj,yj,^j,T]j) , (2.1) 

with 3 real components Xj , yj , Zj and 2 Grassmann variable components , r]j . 
We then define an inner product on K'^'^ by 

(u, u') = -zz' + xx' + yy' + ^77' - 77^' (2.2) 

and constrain Wj by the quadratic equation 

Vj e A : {u,,u,) = -1 , (2.3) 

which is solved by 

z, = ±.^l + x2+y2 + 2^,77, . (2.4) 

In the following we take the positive square root for all j S A . This singles 
out a choice of connected subspace, H^'^, parametrized by two bosonic variables 
Xj , yj and two fermionic variables , r\j . 

On the product space (H^l^)'^' we introduce a 'measure' (more accurately, 
a Bcrczin superintegration form) 

D^i^ = l[{27r)-' dxkdyk d^,d,,, o {I + xl + yl + 2^uVk)-^'^ ■ (2.5) 
feeA 

We use the notation = d/d£, for the partial derivative. The statistical measure 
then is of the Gibbs form Dii\e~^i^-^ with action 

Ap,e ~ - Uj) + £^{zk - 1) (2.6) 

i,j feeA 

= ■^^i (^^^J' ~ (1 + ^i^j + yiVj + ^i^O - Vi^j)) + ^^i^k - ^) ■ 

Here Jij = 1 if i,j are nearest neighbors (NN) and Jij = otherwise. As will 
be discussed later, the action v4^,o is S0(1, 2)-invariant. The e-term breaks this 
noncompact symmetry and makes the integral J Dji\ e~^i^'^ converge. 
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2.2 Horospherical coordinates 

As with , it is very helpful to switch to horospherical coordinates — it is only 
in this coordinate system that we can obtain the probabilistic interpretation of 
Section [L^ We thus use the following parametrization of the supermanifold: 

a; = sinht - e* (ifi^ + V^tA) , y^e*s, ^ c*i/^ , 77 = e*V , (2.7) 

where t and s range over the real numbers. Note that {t, s ^tjjjip) are globally 
defined coordinates and 

(t,s;Vi,V) = (0,0;0,0)^ (a:,2;;C,77) = (0,0;0,0) . 

The expression for the action in them is 

^/'■^ = ^E(.,)(^^. - 1) + ^E.eA^^^ - 1) ' (2.8) 

where (ij) are NN pairs and 

- + - - e*'+*' , (2.9) 

B,j = cosh(t, - tj) + i(si - Sj)2e*'+*^- , (2.10) 
Zk = coshtk + {^sl + i^ki^k) e*" . (2.11) 

We also need the expression for the measure Dfj,\ in horospherical coordinates. 
By applying Berezin's transformation formula [1] for changing variables in a 
(super-)integral, one finds that 

DfiA = Y[^^^{2TT)-^e-*^dtjdsj d^^d^^ . (2.12) 

For any function / of the lattice field variables {tj , Sj , i/jj , 4'j}jeA we now define 
its expectation as 

{f)p,e^ J D^i^e-^^-f, (2.13) 
whenever this integral exists. 

2.3 Effective bosonic field theory 

Since the action (|2.8p is quadratic in the fields ip, ^, and s, each with covari- 
ance Dp^^it)^^ , we know from standard free-field calculus that integration over s 
yields a factor of Det^^^^(£';3_e(<)) while integration over ip yields Det Dp^^{t). 
By performing these integrations, we arrive at the nonlocal free energy func- 
tional Fp^^{t) given by (II. 2|) . Moreover, the basic two-point functions are 

(^PivMv)) = -{[v ; Dp,,{t)-'v]) , (2.14) 
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where 

and the expectations on the left-hand and right-hand side are defined by (|2.13p 
and (|1.5p . respectively. We will often use the formula (|2.14p as well as its 
generalization 

i'iv^mv^)^ ^ _ ^ (2.15) 

where A is the n x n matrix given by 

Axy{t) = K;-D/3,eW"'^^A'] ■ (2.16) 

Remark 2.1. If the Grassmann fields '0,-0 were absent, then Det^^'^ in (|1.2p 
would be replaced by Dct^^^^ (and J^^, e^*''dtk by Hfc e'^''dtk) and Za would be 
the partition function of the hyperbolic sigma model studied in [18] . 

Remark 2.2. If we integrate only over the fields (but not over s) we 
produce a positive integrand depending on t and s. The square root of the 
determinant is then replaced by Det 0/3^^(1) > . 

Remark 2.3. The logarithm of Det_D^.e(t) is convex in t. 

Proof (D. Brydges). By the matrix tree theorem we have 

Det Dp.At) = el^l J] ^^'^'''^ J] . (2-17) 

where !F denotes the spanning rooted forests, R the set of roots, \R\ the cardi- 
nality of this set, and = denotes an edge in the forest. The proof is 
now immediate since any positive sum of exponentials in t is log convex. 

Note that the logarithm of Det Dp^i;{t) competes with the other factor, 
g-/5E(,j)Cosh(t>-tj)^ which is log concave. 

3 Microscopic origin of the model 

In this subsection we use the language and heuristic ideas of physics to sketch 
the origin of our field theory model from a microscopic model of disorder. Con- 
sider real symmetric random band matrices, H, say with independent Gaussian 
distributed entries, of band width W in d dimensions. (Such a band matrix 
model possesses a time-reversal symmetry and belongs to symmetry class Al - 
traditionally referred to as the Wigner-Dyson class of orthogonal symmetry - 
of the 10- way classification of disordered fermion systems TO;). 
Now suppose that we wish to compute the disorder average of 

y/l)et{E + - H)/Bet{E + ie - H) x \ {E + ie - H)-\x,y)\'^ (3.1) 
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for real energy parameters E and £, 77 > 0. The unconventional feature here is 
that the square \{E + ie — H)~^{x,y)\^ of the Green's function is weighted by 
the square root of a ratio of one determinant taken at energy E + ie and another 
one at energy E + ii]. Ahhough one might think that the presence of these extra 
factors comphcates the problem, quite the opposite is true; it will actually lead 
to simplifications when 77 is taken to be large. 
First of all, the combination 

{E + ie- H)-^{x, y) Det^^/^(i; + ie-H) 

can be generated by Gaussian integration over a single real boson field, 0^ . Sec- 
ond, writing the complex conjugate [E + ie — H)~^{x, y) of the Green's function 
as a Gaussian integral requires two real boson fields ^~ and two anticommuting 
fields tjja (Q! = 1; 2). Third, to express the square root of Det(£' + bj ~ H) as a 
Gaussian integral, we need another real boson 0^ and two more anticommuting 
fields Altogether, we then have four bosonic fields (j)'^ and four fermionic 
fields (o- = ±, a = 1,2). 

Now assume for the moment that 77 = e, in which case the two determi- 
nants in p.ip cancel each other. If the band width W is large enough, then the 
standard steps of disorder averaging followed by Hubbard-Stratonovich trans- 
formation and elimination of the massive modes, take us to Efetov's nonlinear 
sigma model for systems with orthogonal symmetry (class Al). 

Physically speaking, the order parameter fields of retarded and advanced 
(— ) type acquire different expectation values: 

(C0^> = 5^p{G^) , (C^?> = ea/3(G^> (a = ± ; a, /3 = 1, 2) , 

where we are using the abbreviations (0„ = (^q^^), 

G+ = {E + ie-Hy^{x,x) , G- , 

and Cap = —e/3a is the antisymmetric tensor for two degrees of freedom. In 
the region of nonzero average density of states, where (G+) 7^ (G^), these 
expectation values break a continuous symmetry of the Gaussian integrand at 
£ = 0. The components of Efetov's sigma model field have the physical meaning 
of being the Goldstone modes associated with this broken symmetry. There are 
4 bosonic Goldstone modes due to the symmetry breaking {(t>a4>t.) {4>Jj4>]}) 
and four more such modes due to {'i't'i't) ~ ~{i'2^t) '^ot being equal to 
{4' I '^2) — ~{'^2i'i)- There also exist 8 fermionic Goldstone modes due to the 
breaking of the odd symmetries connecting the boson-boson sector {<t>1(i>i) = 
(02 "^2 ) with the fcrmion-fermion sector (V'lV'J) — ~{'^2'^i) of opposite type 
T = —a. All these modes organize into a supermanifold with tangent space M^'^ 
over a symmetric space (H^ x H^) x S''. 

Now let 77 ^ e > 0, so that the two determinants in the expression p.ip 
no longer cancel. The difference 77 — e « 77 then acts as a mass term for the 
Goldstone modes connecting the advanced sector (— ) with the 77 retarded sector 
{(t't't't) — '^2) — ~{i'2i'i)- ^ Thouless-type argument, these massive 
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Goldstone modes do not affect the renormalized physics at length scales much 
greater than the length L' determined by the equation 



where D oc is the bare diffusion constant of the system. 

Thus at large length scales L ^ L' we may simply drop the massive Gold- 
stone modes from the theory or, in a more careful treatment, integrate them out 
perturbatively. What we are left with, then, are the 2 + 2 = 4 massless bosonic 
and fermionic Goldstone modes connecting the retarded component ((/>]^0^) of 
the order parameter with its four components {(t)^(t)^) = (^j^^a^) — {ijj^ip2) = 
— (V'^V'r) i'^ the advanced sector. These four residual Goldstone modes orga- 
nize into a supermanifold with tangent space R^'^ and base manifold — we 
thus arrive at the field space H^'^ of the model we are going to study. 

4 Symmetries and their consequences 

As an effective theory derived by reduction from an underlying sigma model, 
the statistical mechanics problem posed by (|l.ip - (|1.5p enjoys a number of sym- 
metries. First among these is a hidden supersymmetry which ensures that the 
partition function is always equal to unity. 



independent of the inverse temperature /3 and regularization parameter e. Thus 
the reduced statistical measure e~^'''^(i/iA can be regarded as a probability mea- 
sure, and the physical observables of the model are given as expectations 



In the following subsection we provide some background to the normalization 
property ZA(/3,e) = 1. 

4.1 Q- symmetry 

We start by observing that, for any e, the full action Ap^^ defined in (|2.6|) is 
invariant under transformations that preserve the short inner product 



for all i,j £ A. Such transformations are given, at the infinitesimal level, by 
even and odd derivations (i.e., first-order differential operators) with the prop- 
erty that they annihilate the expression (|4.ip for all i,j and their coefficients 
are linear functions of the coordinates Xk,yk,£.k,Tlk ■ These differential opera- 
tors form a representation of the orthosymplectic Lie superalgebra osp2|2 • 
important example of an odd operator Q e osp2|2 



r/ = 2-KhD/L' 



ZA{(5,e) = 1 




(4.1) 



(4.2) 
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Since Y[j dxjdyj d^-dri^ is the Berezin superintegration form given by the inner 
product (j4.ip . it is immediate that Dn\ is osp2|2"ii^^9'^i8''^t, which impHes that 
/ D/iA Qf = whenever the function / is differentiable and Qf is integrable. 

For present use, let us record here the explicit expression for the osp2|2 S^n- 
erator Q in horospherical coordinates: a straightforward computation starting 
from (|4.2p gives Q = J2j£A9j with single-site generator (index j omitted) 

q = ^dt + ii^- s^P)ds -sd^ + ^{l- e-2* - s2 _ i^^)d^ . (4.3) 

Now consider any differentiable integrand / which is invariant by Q, i.e., 
Qf = . This invariance property has strong consequences for the integral of / 
(if it exists): in Appendix O Proposition [51 we prove that the integral of such / 
equals / evaluated on the zero-field configuration (i.e., on tj = Sj = -tjjj — -tpj — 
or equivalently, Xj = yj = — Vj = ; for all j G A) : 

JD^lAf^fio). (4.4) 

The idea of the proof is easy to state: one shows that the integral of / remains 
unchanged by the replacement / c^^'^f with h = J^jeA^^^ + + '^^jVj) ^^'^ 
T > , and then deduces the result (|4.4p by sending the deformation parameter 
T — » -|-oo to localize the integral at the zero-field configuration. 

Using the explicit expression (|4.3[) it is easy to check that the action Ap^^ 
is Q-invariant. Since the differential operator Q is of first order, one directly 
infers the relation Qe""^"-^ — 0. Therefore, as a particular consequence of ()4.4|) 
and Ai3^^{o) = it follows that the partition function equals unity, 

Za(/3, e)= I D^iA e-^^- = e'^^-^") = 1 , (4.5) 

for all values of /? > and e > 0. 

Further consequences of (|4.4p will be elaborated below. 

4.2 Hyperbolic symmetry 

While Q is a symmetry of our action Ap^^ for all values of e, further symmetries 
emerge in the limit of vanishing regularization e 0+. Relegating a more 
detailed discussion to Appendix [Bl we here gather the crucial facts. 

The model (|l.ip - ()1.5p for e — > 0-1- acquires a global symmetry by the Lorentz 
group S0(l,2) - the isometry group of the hyperbolic plane viewed as a 
noncompact symmetric space ~ S0(1, 2)/SO(2). This global symmetry 
entails a number of conserved currents and associated Ward identities. Of these 
let us mention here the most important one, 

J2{e'^+'yDp,,it)-\x,y))^-^, (4.6) 
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which is the sigma model version of the quantum sum rule 

J2{\iE-ie- H)-\x, y)\^) = \ (Im (£; - ie - H)-^{x, x)) = -^p{E) , 

where p{E) is the mean local density of states. In the sigma model approx- 
imation one sets Trp{E) = 1. The above relation reflects the unitarity of the 
quantum theory. Its classical interpretation is conservation of probability. 

Notice that the right-hand side of (|4.6p diverges in the limit of vanishing 
regularization e — > . For an infinite lattice A there exist two principal scenarios 
|12| by which to realize this divergence. In the first one, the correlation function 
Cxy = (e*^+*!'£>^^(i)(x, ?/)), while bounded in e, becomes of long range and 
thus ceases to be summable in the limit e ^ 0. In this case the S0(l,2) 
symmetry is spontaneously broken and the system is in a phase of extended 
states. On the other hand, Cxy may already diverge for any fixed pair of lattice 
sites x,y, signaling strong field fluctuations and restoration of the noncompact 
symmetry S0(1, 2) as e — > . Exponential decay of Cxy with distance \x — y\ 
then corresponds to exponential localization of the energy eigenstates. Thus the 
question of extended versus localized states of the disordered quantum system 
translates to the question of the Lorentzian symmetry S0(1, 2) of the statistical 
mechanical model with free energy being spontaneously broken or not. 

At this stage, a remark is called for. Niedermaier and Seller have recently 
shown [m [16] for a large class of sigma models that if the symmetry group 
of the sigma model is non-amenable - this includes in particular the case of 
the Lorentz group S0(l,2) - then spontaneous symmetry breaking occurs in 
all dimensions d > 1 and for all /? > . It must therefore be emphasized that, 
although our sigma model does acquire the non-amenable symmetry S0(1, 2) in 
the limit e ^ 0, it does not belong to the class of models where the arguments 
of [m [16] apply. (The culprit is the nonlocal part of the free energy due to 
integration over the Grassmann fields.) In fact, the non- amenable symmetry 
S0(1, 2) of our model is known to be unbroken in d 1. This follows from the 
work of [21] where the conductance of the one-dimensional system was shown 
to exhibit exponential decay with increasing length of the system. 

4.3 Perturbative renormalization group 

We now sketch a perturbative result from Wilsonian renormalization theory 
by which our model is expected to be in a symmetry-unbroken phase also for 
d ~ 2 and all values of the inverse temperature /?, and thus to exhibit Anderson 
localization of all electronic states. 

This result follows from Friedan's work [8] on renormalization for the general 
class of nonlinear sigma models. According to it, the RG flow of the temperature 
T = (3~^ with increasing renormalization scale a is given by 

= [2 - d)T + RT^ + 0{T^) , (4.7) 

da 



13 



where R is the target space curvature - more precisely, the multiphcative con- 
stant R by which the Ricci tensor of the target space differs from its metric 
tensor. For both the H^'^ model and Efetov's sigma model of class Al a quick 
computation shows the curvature R to be positive. In contrast, R = and 
R < for Efetov's sigma models of class A (broken time-reversal symmetry) 
and class All (spin-orbit scattering), respectively. 

According to (|4.7p . a positive value of R implies that a small initial value 
of the temperature T increases under renormalization in dimension d = 2. By 
extrapolation, one therefore expects the existence of a mass gap (or, equivalently, 
localization of all states) in this case. For the localization length ^ = ^(a, T(o)), 
which is a physical observable and hence a renormalization group invariant, one 
obtains the formula 

Ccxae^/^^^) (d = 2) 

by direct integration of the RG equation (|4.7p . 

In dimension d = 3, equation (j4.7p predicts the localizing tendency of posi- 
tive target space curvature to become irrelevant at small enough temperatures 
and hence the RG flow to be attracted to the fixed point T = corresponding 
to extended states. As was remarked above, the Lorentzian symmetry S0(1, 2) 
is spontaneously broken at this fixed point. 

With increasing temperature T (or decreasing field stiffness /?) the model 
in d = 3 is expected to undergo an Anderson-type transition to the phase of 
unbroken symmetry. This phase transition was studied numerically in [4] , where 
the critical value of /3 was found to be (3c ~ 0.04. The transition has also been 
investigated in detail using the Migdal-Kadanoff renormalization scheme [5] . 

5 Ward identities and outline of proof 

In order to control fluctuations of the field t at low temperatures T = (3^^ 
we rely on a family of Ward identities due to the internal supersymmetries of 
the model. These Ward identities are naturally expressed in terms of both the 
real variables tj , Sj and the Grassmann variables ipj , ipj . In order to obtain 
probabilistic information we integrate out the Grassmann variables using (|2.14p 
and (|2.15p . thereby producing a Green's function. 

As was already mentioned, our partition function always equals unity even 
when the temperature varies in space. By using this fact, we show that gradi- 
ents of the field t between neighboring sites are strongly suppressed for small 
T. There also exist Ward identities at larger scales, and information may be 
extracted from them by using information on previous length scales. 

In addition to Ward identities, there are two other crucial ingredients of our 
proof. The first one is a basic estimate on Green's functions which are non- 
uniformly elliptic. The second one is the use of SUSY characteristic functions, 
which help to control large-scale field fluctuations. 

A more detailed outline of our proof is given below, where the notation 
and the needed Ward identities are explained. Once the Ward identities are 
established, most of our proof is very classical. 
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5.1 Ward identities due to Q-symmetry 

We recall the formula (|4.4p for the integral of a Q-invariant function. It is 
easy to check that the functions Sij and given in (|2.9p and (|2.1ip satisfy the 
invariance conditions Q5y = and Qzk = 0. Therefore, using Sij{o) — 1 and 
Zk{o) = 1 we have the identity 

J DfiA e-''^-.«"'==''(^==''"^)"S-eA^-(^--i) = 1 (5.1) 

for all values of > 0, Jxy > 0, and Sx > 0. Note that in order for this statement 
to be true, J^y does not have to be nearest neighbor. Moreover, for m £ M and 
any pair x,y d A we have 

i = (^."W = (s;';+"^^r;-'(^.-^«)(^.-V',)e*-+''')^^^ , (5.2) 

where the expectation {■)(3,e was defined in (|2.13[) . and we used the nilpotency 
{ipx — ^yY = {i'x — i^yY — 0- By integrating over the Grassmann fields i/; and 
ip as in (|2.14p we obtain our basic identity, 

l = (^r>e = (i?:;(l-mG,,)). (5.3) 

The last expectation is taken with respect to the effective action for the fields t 
and s, and the Green's function Gxy is 

Gxy = [{.5x - Sy); Dp^Sr\5x - <5,)]^ . (5.4) 

iJxy 

More generally if {xi,yi) are n pairs of points, then 

1 = (IIL ^"..>,,, = KyP-^i^ - -^)) > (5-5) 

where Q is an n x n matrix of Green's functions 

g^J^[g^;Df3Atr'gJ]A (5.6) 

and 

S.^S.-.lf e(*-+*-)/2(<5..-<5,J. (5.7) 

The matrix Q is real symmetric and positive. It will be important later that we 
are choosing gi to be orthogonal to the zero mode (i.e., the constant functions). 

5.2 Outline of proof 

Our proof of Theorem [T] relies on the Ward identity (|5.3[) . (|5.5p and an induction 
on length scales. The basic idea is quite simple: suppose m > and we had 
a uniform bound \Gxy\ < C//3 < 1/m on the Green's function ()5.4|) . for all 
configurations of t. Then we could conclude from (|5.3p that 

(cosh™(t, - ty)) < {B^y) < (1 - mC//3)-i , (5.8) 
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Figure 1: "diamond" region: a double cone in 3 dimensions 

and this would imply Theorem [TJ 

In Section [S] we prove that if — y| = 1 then indeed < Gxy < i/f3 , and we 
establish an even stronger version of (|5.8p . This proves that nearest neighbor 
fluctuations of the field t are very unlikely for large /3 (see Lemma |4]). 

For distances |x — j/| > 1, however, there is no uniform bound on Gxy ■ In 
Section [7] we study the Green's function (|5.4p and estabhsh sufficient conditions 
on the field t to obtain the desired bound on Gxy ■ In 3D these conditions are 
roughly given as follows (where |j — x| > 1): 

cosh{t J -tx) < Bjx < a\j - x\" , 0<a<l/2, (5.9) 

and the same for cosh(tj — ty). The number a is a constant, say a > 10. It will 
turn out that these estimates are needed only for the sites j in a 3D diamond- 
type region, Rxy , containing x and y ; see Fig. [T] Notice that since the exponent 
a is positive, we are allowing larger fluctuations at larger scales. The probability 
that such a condition is violated will be shown to be small by induction. 

Section [8] uses the conditions described above to prove conditional estimates 
on the fluctuations of the field t at all scales. These conditions are initially 
expressed in terms of Q-invariant characteristic functions x- Later we show 
that the nilpotent (or Grassmann) part of x is not important, so we may think 
of X in the usual classical sense. 

The remaining problem is to obtain unconditional estimates on the fluctua- 
tions and thereby prove Theorem [T] This is first done for short scales in Section 
[TOl For larger scales we use induction. Our induction hypothesis is formulated 
in Section [TT] Roughly speaking it asserts that 

(nL^":..>^2- (5.10) 

holds under the assumption that the diamond-type regions Rxiyi (Fig. [l]) asso- 
ciated with i = 1, . . . ,n have disjoint interiors. The induction is in £, defined as 
the maximal separation |a:i — in the product over i — 1, . . . , n. For £ ^ I this 
hypothesis was verified in Section [6l 
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Section [T^ contains the technical core of our paper. There we prove uncon- 
ditional estimates on the fluctuations and thus obtain Theorem [TJ The main 
idea is to consider a site b in R^y closest to a; or y such that condition (|5.9p is 
violated for j = b. We shall then prove by induction that the probability for 
such an event to occur is small. The inequality i3™ < 2'"_B™B™ (see Lemma 
[5] below) is used for a point c near b. Since the distances \x — c\ and |c — y\ are 
less than \x — y\, induction can be applied. The factor 2™ is offset by the small 
probability of the event when /3 is large. 

Theorem [5] is proved in Section [T51 Here we must estimate the contribution 
of the zero mode and at this stage e > plays a key role. Finally, Theorem [3] 
follows from the estimates of Theorem [2] its proof is given in Section [TH 

5.3 Two simple lemmas 

We conclude this section with two simple lemmas which will be frequently used 
below. The first lemma is useful for estimates on Green's functions. To state it, 
let y be a finite-dimensional Euclidean vector space with scalar product [ ; ]y. 

Lemma 1 Let M : V ^ V be a positive real symmetric operator. Then for 
any set of n vectors Vi €z V we have 

n 

M~^t;,K; -Iv >0 (5.11) 

'i=i 

if and only if the symmetric n x n matrix K with matrix elements 

K,j^[v,-M-\,]v (5.12) 

satisfies < X < Id. 

Proof. Letting Wi = M^^^'^Vi € y we observe that Kij — [wi ; 'Wj]v and (|5.1ip 
is equivalent to 

n 

ld-^^Wi[wi; ■]v >0 ■ 

1=1 

By evaluating this quadratic form aA, w — ^ for any real numbers Xi we 
see that (|5.1ip is equivalent to 

n 

J2 >^^>^J {K ~ K%, > , (5.13) 

or < < K, from which our assertion follows. ■ 
The second lemma will be used in our induction process of Section 1121 

Lemma 2 // B^y , Bey , B^c o,fe defined as in (j2.10p then 

Bxy < 2BxcBcy ■ (5-14) 

The inequality ([STM]) can be verified by direct computation (proof omitted). 
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Remark. The raison d'etre behind (|5.14p is easy to state: B^y has an inter- 
pretation as the hyperbohc cosine of the geodesic distance on H^. Therefore, if 
x,y,c are three points on H^, then since the geodesic distance dist(x,j/) is the 
minimal length of any curve connecting x and y, the triple of geodesic distances 
satisfy the triangle inequality dist(a;, y) < dist(a;, c) + dist(c, y). Given this, the 
inequality (|5.14p follows by taking the hyperbolic cosine of both sides and using 
that cosh(a + b) < 2 cosh(a) cosh(&) holds for any two real numbers a, b. 

6 Bounds on NN fluctuations 

As was already mentioned, for nearest neighbor (NN) pairs we can obtain a 
result stronger than (|5.8p . Recall that we have now fixed Jxy = 1 for all xy that 
are NN pairs, and Jxy = otherwise. This fact is essential in the next lemma. 

Lemma 3 Let x, y he an NN pair and suppose that < 7 < 1. Then 

(e'^''(^--i))<(l-7)-i. (6.1) 
More generally, if {xj,yj), j — 1, . . . ,n is a set of n different NN pairs, then 

^e/37E?=i(Sx,-.,-i)^ < (1-7)-" . (6.2) 
This shows that NN fluctuations are strongly suppressed. 

Remark. Since Jxjyj = 1 and 7 < 1 the integrals in (|6.ip - (|6.2[) are well 
defined. This would not be true if Xj, yj were not NN, or if 7 > 1, or if two or 
more NN pairs were allowed to be identical without further restrictions on the 
value of 7. 

Proof. For x, y an NN pair let 

Fxyij)^e^'^+'y^/'i6xij)^Syij)), (6.3) 
and introduce the Green's function 

Glyit) = [Fxy ; D^^^it) Fxy] = BxyGxy . (6.4) 

Since Sxy is Q-invariant, Proposition [2] of Appendix [Cl implies 

g/37 = ^g/37Sxy^ ^ ^g/37(Bxy+^'(Fx„)^(i=^xy))^ 
^^^^P^B^y^l+P^^i^FxyMFxy))) 

= (g/37Sxy(i_^^GS^))^ (6.5) 
where we have used i/P = = and (|2.14p . Now from (|l.ip we have 

[v ; D0^,{t) v] = (3j2u ■) ' ^^^1' + ^ ^ ^ ; ^-y^' ■ (6.6) 
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Therefore Lerama [T] implies that < f3G'^y(t) < 1 for all t, and (|6.ip follows. 
Similarly, for 71 > 1 we have 

= ^e/37E^=iB.,«,Det(l-7if)^, (6.7) 

where K is the n x n matrix 

K^, = /3 [F.,^y^ ; i^/3,e (i)-' F.,„J (6.8) 

given by n different NN pairs Xi,yi . From (j6.6p and Lemma [T] it follows that 
\\K\\ < 1. This implies |Det(l — ^K)\ > (1 — 7)" and the lemma follows. ■ 
As a corollary, since 1 < < e™(^-«-i) for m > 0, we have the bound 



n,=i^"..)^(l-"^//3)""<2" (6.9) 

for any m in the range m < /3/2 . 

A first important consequence of Lemma [3] is the following statement. 

Lemma 4 Let Xj, yj be a set of n different nearest neighbor pairs. Then 

Prob (Vj = 1, . . . , n : B^.y^ > 1 + ^) < (1 - 7)-"e-"('37)5 (g.io) 
for any < 7 < 1. 

Proof. Let n = 1. By the Chebyshev inequality , 

Prob(B,, > l + <5) = (x(B,, > l + ,5)) <e-'"^(i+^)(e''^^-) , (6.11) 

where xi^xy > 1 + 5) is the characteristic function for B^y > 1 + S to hold. The 
desired inequality for n = I now follows directly from Lemma [3] 

The proof for n pairs is no different. ■ 



7 Conditional estimates on Green's functions 

For general x,y (not NN) we do not have the option of considering (e^''''^^«), as 
the underlying integral need not exist. Nevertheless, (5*^) does exist and from 
(|5.3p we have 

l^{Bll{l-mGxy)) , (7.1) 

with Gxy defined by 

Gxy = [{Sx - 5y);Dp,e{t)~'^{5x ~ Sy)] . 

^Actually, the Chebyshev inequaUty states that for any random variable X with average 
Jfo we have Prob[(X — Xq)^ > a^] < ((X Xq)^). Here we are using the same principle. 



19 



Now, as was explained in Section 15.21 if we knew that Gxy ^ C/(3 for all con- 
figurations of t, then we could conclude that 

TO<(l-™C//3)-i. (7.2) 

While we have seen that this estimate is true for |a; — y| — 1 (with C = 1), it is 
false in general, as there are rare configurations with large negative t surrounding 
X or y. Nonetheless, in 3D we can get an upper bound on Gxy by estimating 
the local 'conductance' at an edge {ij) from below. This conductance is 

> imax(e*'+*^-2*^,e*-+*^-2*«) , (7.3) 

where we have used Bxy > cosh(t2, — tj,). It will sufRce to estimate the expression 
(|7.3p for NN pairs (ij) in a region Rxy which is like a 3D double cone with vertices 
at X and y. Note that Neumann boundary conditions increase Gxy and Ox — Oy 
is orthogonal to the zero mode. We will have to require that Rxy be essentially 
three-dimensional in the following sense: 

Definition 1 A region Rxy C A containing x and y is called S-admissible if it 
is connected by nearest-neighbor bonds and the two one-parameter families of 
intersections Rz{r) = Rxy H B^. with the ball B], of radius r centered at z = x,y 
satisfy 

vol{Rz{r)} > r^6 for r<\x — y\/V2 {z = x,y). 
In addition we require that the following Poincare inequality: 

E /(jf< const (W)'W' 

holds for all functions / : A — s- R subject to the condition X^jei?. (r) f(j) ~ ^ ■ 

We observe that by the choice of maximal radius r ^ \x — y\/V^ the scaling of 
volume is monitored up to the full side length of a rectangular diamond Rxy (or 
a double cone Rxy , see Fig. [TJ with opposite corners placed at x and y. 

In the continuum limit this definition is satisfied by a double cone obtained 
by rotating (around the line xy connecting x and y) a 2D diamond with vertices 
on X and y and angle 9 > 9o{S) > tt/IO (see Fig. [5^). Since we are on a lattice 
we may have to add a few lattice points near x and y to ensure connectedness 
(see Fig.lJla). The Poincare inequality is straightforward to prove in such convex 
regions. 

Definition 2 Given a S-admissible region Rxy , we define the regions Rxy for 
z = X and z = y by 

Rly = {3 e Rxy \l<\3~z\<\x- y\/V2} . (7.4) 

For the case of a diamond, R!^ U R^. — Rxy \ {x, y} . 
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Figure 2: a) The region Rxy in the continuum limit, b) On the lattice, the 
points xi and X2 must be added to ensure connectedness. 

Remark 7.1. The values of the field t outside the region R^y are not impor- 
tant, as we can use Neumann boundary conditions to eliminate the exterior of 
Rxy ■ Indeed, in the subspace orthogonal to the constant functions the Lapla- 
cian on R^y with Neumann boundary conditions is bounded (by the Poincare 
inequality) from below by some number, say c , times the inverse square of the 
linear size L of Rxy ■ By this token, since the vector 6x — 5y used in the defi- 
nition of Gxy lies in that subspace, we may utilize the bound on the inverse of 
the Neumann Laplacian by c~^L^ and in this way eventually obtain an upper 
bound on Gxy (see Lemma [5]). 

Remark 7.2. From (|7.3p we have 

^xyiij)'^ < 8 cosh(ti — tz) cosh{tj — tz) (7.5) 

for both z ~ X and z = y. The main result of this section is that we can 
get an upper bound on Gxy even without imposing a uniform upper bound on 
Axy{ij)~^, as long as some growth restriction on the fluctuations of t is met for 
a (5-admissible region Rxy : 

Lemma 5 Fix two constants a> 1 and 1/2 > a > . If Rxy is a 5-admissible 
region in the sense of Def. Q] and the statement 

Vj e Rxy ■■ cosh(tj -tz) < a \j - z|" (7.6) 

holds for both z — x and z — y, then we have 

0<Gxy<G^y<C{a,a,6)/p , (7.7) 

where C{a, a, S) is some constant depending on the parameters a, a and the 
geometry of the region Rxy (encoded in the parameter 6). The notation G^ 
means Neumann boundary conditions on dRxy ■ 
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Proof. The following is a variation on an argument presented in [IS] . For each 
A; G N consider two cubes of side 2^ centered at x and y. (For concreteness, 
imagine the edges of the two cubes to be parallel to the vector x — y.) Let 
, Ry denote the corresponding intersections with R^y and let Ik , Ik be the 
indicator functions of R^ and Ry , respectively, normalized so that for each k 

5].4(j) = l = 5]./fe(j). (7.8) 

We observe that 

yk > kra : i?^ = R^ = Rccy , (7.9) 

where km is the smallest number fc G N such that 2*^^^ > |a; — yj. Since R^y 
is (5-admissible, i?^ has the same properties as Rz{r — 2*^) in Dcf. [1] and we 
therefore have 

voli?^ > 2'^'*(5 {k<k„,-l, d = 3) (7.10) 

for all 2'^ > 10. For 2*^ < 10 this is not true (see Dcf. [T]) but the corresponding 
volume is no less than unity, as R^y is connected. 
Now we express Sx — 5j, as a telescopic sum: 

S^-5y = ^ {pk - Pk) , (7.11) 
fc=i 

where pk = h-i - h , Pk = h-i - h , h = 5x , and Iq = 5y . This sum 
terminates at km because by (|7.9p we have Ik = Ik for k > km ■ Note that pk , 
pk are orthogonal to the constant functions: X^jeA PkiJ) = = '^jeA PkU)- 
Next we put the telescopic sum to use by the following computation: 

km 

[Sx-Sy;Dp^,{tr\5x-5y)] = J2 [pk-pk;Df,Atr\pi -Pi)] 

k,l=l 

^ (E(^fc;^/3-^(*)"V/c]'/'+[pfc;i?/3,e(i)-^P.]'^')^ , (7.12) 

where the Cauchy-Schwarz inequality was employed. Hence we need to estimate 
[pk ; D/3^e{t)^^ Pk] and [pk ; Df3^e{t)~^ Pk]- This is done, say for the former, by the 
inequality 

[pk;DpAtr'pk]<\\D^l\\\\pk\\l, (7.13) 

where Dfi (for a region R) stands for the operator (jl.3p with Neumann boundary 
conditions on R. In view of ^ Pk{j) = the operator norm is to be taken on 
the orthogonal complement of the constant functions. 

The square of the L^-norm of pk+i is bounded by (voli?^)^^. Thus by (|7.10p 

\\pk+i\\l<{yolR^)-^<2-''''d-' . (7.14) 

The corresponding inequality also holds for pk+i ■ 
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We must still bound the operator norm \\D„^. \\. For this we observe that the 
conditions (|7.5p ensure that 

for all j, j' e and k<km — l, since in that case jj' G R^y and we apply (|7.6p 
for z = X. For = km we are looking at pairs that belong to i?^^ but not to 
R^y. In that case we apply (|7.6p for z = y and still have A^yijj')^^ < 80^2^*^". 
Therefore, since Rxy is (5-admissible and the lowest nonzero eigenvalue of the 
Neumann Laplacian on R^ is of the order of (2*^)"^, we obtain 

e'^+*^B-^\\D-^l\\ <c{5) 0^(3-^2^'' 2^'"^ (7.15) 

for some c{5) and all k < km ■ For 2*^ < 10 the connectedness of Rxy ensures 
that \\D~l\\ < const. The same bounds apply for ||Dd;,||. 

Finally, by combining ((77121) with (fTTSl) . jZHl]), and (fTTSl . we arrive at 
G.. = ^ ['J. - ; Dp,,(t)-\5x ~ 5y)] < c{6) ^ ('2£ V2'^'"+'^'') . 

'■^ y ^ \ k=i J 

For 2a < d — 2 = 1 the value of this sum is bounded imiformly in fcm . ■ 



Remark. The bound (|7.7p also applies when the definition of R^y and is 
modified in the following way (for z — x,y a,s before): 

- {j e Rxy : |j -z\<\x-y\f,} , (7.16) 

where fx , /y are a pair of positive numbers which add up to (at least) unity and 
neither of which is too small. It is easy to see that the relevant scales involved 
are the ones for k near fc„i and we can get the same bound but with a change 
of overall factor. This remark will become important in Section [121 Lemma fTTl 
where we will need this estimate with /j, ~ 1/5 . 



8 Conditional estimates on fluctuations 

In this section we establish bounds on the fluctuations of the field t by bounding 
(i?™ Xxy) where Xxy has the property that Xxy = whenever the hypothesis 
(|7.6p of Lemma |S] fails. 

Definition 3 (characteristic function); As before, fix two constants a > 1 
and 1/2 > a > 0, and let rj-k '■— (a |j — fc|")~^ for j, k £ K, j ^ k. Let 
X '■ — > M fee the characteristic function of the interval [0, 1], i.e., x{t) — 1 
for < ^ < 1 cind x(i) = for t > I. Moreover, let Rxy be S-admissible and 
choose the regions Rxy , R%y as in (|7.4|) . In this setting we define 

Xxy= n n ^ai ' Xzj =x{'r]-zB^,) iz^x,y). (8.1) 
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Here the constants a, a are taken to coincide with those in Lemma 
With these definitions we have 

Lemma 6 Let R^y be a 5 -admissible region, and let C = C{a,a,S) be the 
constant that appears in Lemma [31 Then for < m < P/C we have 

{B^yX.y)<{^-mC/f3r' . (8.2) 

Proof. Our proof uses the identity (S*™ Xxy) = 1 where xty is a supersym- 
metric version of Xxy defined above. After integrating out the Grassmann fields 
we shall show that this identity implies 

{B-^yX.y{l-^Gxy))<l. 

Lemma [S] and the presence of Xxy then yield (|8.2p . 

More precisely, let Xi € C°°(M+) with j-^x-yit) < and 



X7W = 




t< 1-7, 
t> 1 , 



be a smooth regularization of x = lim^^o X-y ■ We fix a small value of 7 > and 
write x = x-y for short. Then, recalling the definition (|6.3|) of Fxy we introduce 

Xxj = Xirj-x Sxj) = Xxj + rj^x Xxj ^iFxj)^iFxj) , (8.3) 

where Xxj — xi'^'j-xBxj)- Since Xxj Q-invariant and xfj(O) = 1 we have 

1 = (Sxy n n 4) ^ (STyX'xy) ■ (8-4) 

Now, we express 

{SZ Xxy) = {STy Xxy exp - ; AV-]), (8.5) 
where the symmetric operator A is given by 

[f;Af] = - ^ 'JZ^[f-Fx,]' - E '-^If-^FyA' ■ (8-6) 

Clearly ^4 > as a quadratic form since x' < . The total contribution to 
(|8.5p including the fermionic part of the action is 

; [Dp^.it) + A)^] - mB^y'^PiFxy) i^{Fxy) , 

where the second summand stems from 5*™ , see (I5.2l) - (|5.3p . Thus, integration 
over the Grassmann fields ^/i, i> gives T)ei{Q + A) where 

Q = Dp,,{t)~7nB-^Fxy [Fxy; •]. 
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Since we are taking m to be less than /3/C, the presence of the factor Xxy in 
(|8.5p ensures (by Lemma [5]) that mGxy < 1 • Now by Lemma [T] the inequahty 
1 > mGxy = 'mB~y[Fxy;Dp,^{t)~^Fxy] is equivalent to Q > . Therefore the 
result Det{Q + A) of integrating over tp, is bounded from below by 

Det(Q + A)> Dct(O) = Det(D^,£(t)) (1 - mGxy) > , 

and we obtain the estimate 

1 = {BTy X.y e-t^ > (S™ x.y (1 - niGxy)) • 

We finally take the limit 7^0. The smooth function Xxy then converges to 
the characteristic function Xxy ■ Hence 

1 > (Sr, Xxy (1 - mGxy)) > {B^:y Xxy) (1 - mG/P) , 
which is the desired result. ■ 



Lemma 7 // all of the regions Rxiyi 1 Rx2y2 ; ■ ■ • ^ B,x„y„ are 6-admissible and 
disjoint (meaning they have disjoint interiors), then we still have 



n 



B:^^y^Xx,y,) <{l-mG/P)-- . (8.7) 



Proof. As before we use the fact that the supersymmetrized observable, which 
here results from replacing Bx^yj by Sx^y^, has expectation one. 

Consider first the simpler problem of computing the expectation of the prod- 
uct Y\S^xiyiXxiyi ■ After integrating over ip and ^ we see that 

^n^^y, xx.v}^ = ^n^"?.. ^-^y^ - '"^)^ ' (8.8) 

where ^ is an n x n matrix of Green's functions 

g,, = /3-i[g. ; Dp^,{t)-^g,] , g, = i?-^^ e(*^.+*='.)/2 (Sx, - SyJ . (8.9) 

The matrix Q is positive as a quadratic form. In order to reduce the problem to 
the previous case (of just a single region) note that G < Gn , where the subscript 
denotes Neumann boundary conditions on the boundaries of the disjoint regions 
Rxiyi ■ The presence of the factors Xxiyi implies bounds on the Gx^y^ so that 

n 

Det(l - mG) > Det(l - mGw) ^ ^[{1 - mGx^y,) > (1 - mC//3)"" . (8.10) 

i=l 

The proof of the lemma is completed by introducing the effects of x' 
before. Since there are no new aspects to this argument, we omit it. 
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Remark. From this lemma one obtains estimates for conditional probabilities 
only. Yet, in order to bound Cxy in Theorem [3] we need probability estimates 
without any conditions, which is why we now have to develop an inductive 
argument. 

9 Unconditional estimates on fluctuations 

We are now going to remove the constraints enforced by insertion of x- In order 
to do so, we have to consider x%j — ^ ~ Xxj for Xxj defined by ()8.ip . Short 
scales (given by < \ j — x\ < (3-^^^) will be treated separately by monitoring, in 
Section [ini only the size of nearest neighbor gradients inside the region Rxy • -^t 
the very large scales of \ j — x\> /3^/^, however, looking only at NN fluctuations 
is not enough. There, in order to remove the x constraints we will show by 
induction on the distance \j — x\ that the corresponding contribution is small. 

We will distinguish between two types of geometry: diamonds and deformed 
diamonds. For deformed diamonds we will quantify the bounds given by (|7.7p 
and call such regions C-admissible. 

Definition 4 Let R^y d A be 5 -admissible in the sense of Def. [H 

1. We call Rxy a diamond if it is the set of lattice points which is contained 
in a 3- dimensional double cone obtained in the following way: we take a 
2-dimensional rectangular diamond with opposite vertices placed on x and 
y and edges of length \x — y\/^/2 (see Fig.\B^ and rotate it around the line 
xy. In order to ensure connectedness we may have to add a few lattice 
points near x and y (see Fig. \^a,b). 

2. We call R^y a C-admissible region (or deformed diamond) if 

0<G^yX.y<C/l3 for\x-y\>(3'/\ (9.1) 
0<G%j]Xj,g^P^^^/f^ /or - 2/1 < (9.2) 

where Xxy is defined in (|8.ip . the superscript N stands for Neumann 
boundary conditions on R^y and J^' denotes the product over all near- 
est neighbor pairs in Rxy ■ 

Note that for short scales, dealt with in (|9.2p . instead of using Xxy we impose 
constraints on all NN pairs in the region Rxy ■ 

With these definitions we can state the main result of this paper. 

Theorem 4 Let m = (3^^^ , and let Rxiyt for i = 1, . . . ,ni be diamonds with 
disjoint interiors and \xi — yi\> (i^l^ . Then we have 
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a) b) 

Figure 3: a) diamond region, b) C-admissible region (deformed diamond) 

for all Til > 0. Moreover ifpj, qj for j — 1, . . . ,n2 are such that \pj — qj \ > /?^^^, 
the regions Rp^q^ are C-admissible, have disjoint interiors and do not overlap 
with any of the regions Rxiyi , then there exists a constant p < 1/2 such that 

n^"..n^m.xp.^.) ^2-(i+p)"^ (9.4) 
i=i ]=i I 

for all ni > and 712 > 0. Finally, let rk,Sk for k — 1, ... ,713 be such that 
\rk — Sk\ < P^^'^, RrkSk o-f^ C-admissible, have disjoint interiors and do not 
overlap with any of the regions Rxiyi or Rp^q^ ■ Then for all rti > 0, ri2 > and 
> we have 

(rii n2 "3 \ 

n n n ^rix ) < (i + pr 2- (9.5) 
1=1 J"=l fe=l / 

with p < 1/2. 

The proof of the theorem is carried out in Sections [T0l - [T2l We will need to 
distinguish between three situations, which we refer to as classes. 

Class 1. \x — y\ > P^^"^ and the pair is not protected by a factor of Xxy ■ In 
this case we need an inductive argument on scales to prove a bound on the 
expectation of i3™ . The induction will be done on £ — max^ \xj — yj\ and 
is carried out in Sections [TT] and [T^ We will need to inductively select non- 
overlapping smaller diamonds inside the region R^y while making sure that 
these remain (5-admissible. To arrange for all geometrical details to work out, 
we take Rxy to be a perfect diamond. 

Class 2. |a; — ?/| > P^^^ but the pair is protected by a factor of Xxy ■ In this 
case we can apply the results of Section [8l thereby obviating the need for any 
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induction. R^y is then allowed to be a deformed diamond and the bound we 
can get is stronger than in Class 1 (power 3m instead of m). 

Class 3. \x — y\ < f]^/"^ . This includes short scales and the NN case, which was 
already treated in Section [S] We will show in Section [TO] that these scales do 
not require any factor of Xxy to ensure a good bound. No induction is needed, 
and we can therefore take R^y to be a deformed diamond. 

Note that the larger exponent 3m appearing in (|9.4p and (|9.5[) is important 
for the inductive proof to go through. The enlarged exponent can be handled 
either because of the presence of x or because the pair is of Class 3. 

9.1 Fixing the different parameters 

We have introduced a certain number of parameters: m, a, p, S, C , a. Before 
going on, we briefly review why they appeared and how to choose their values. 

1. The parameter m is ubiquitous in this paper as the power of Bxy ■ Since 
the probability of large deviations will be bounded by K^™ with K > 1, 
we want m to be as large as possible. On the other hand, to apply the 
SUSY argument of Section [8] we must have (1 — 3mC//3)~^ < 1, where 
the factor 3m in this inequality comes from the power of B in (|9.4p . (|9.5p . 
Therefore the magnitude of m is limited by (3. To arrange for all the 
conditions to be met, we fix m = 0^^^ . The factors m and 3m will be kept 
fixed in the whole course of proof. 

2. The constants C and 5 appearing in the definition of the region Rxy (see 
Def. [T]) are not subject to any special requirements, but their values do 
constrain the other parameters. They will be fixed throughout. 

3. To prove the induction hypothesis we need < p < 1. More precisely (see 
Eq. ()12.6p ) we need p + Tl{x, y) < 1. Since we prove TZ{x, y) < 1/2 we will 
take p < 1/2. 

4. The constant a in Lemma [5] in Section [7] plays a key role in bounding 
the entropy for small scales; see Section [T^l Case 1, Eq. (|12.16p . It will 
become clear there that a > 10 is sufficient. 

5. We need to take a > in Lemma [5] in order to control entropy factors 
for large deviations (see Section [TO)) . On the other hand, the result of 
Theorem 2 would be optimal for a = 0. Therefore we wish to make a as 
small as possible. We will see in Section [TO] (Case 2b eq. (|12.2ip and Case 
2c eq. (|12.23p ) that a > 0(1/ In /3) is a requirement for our analysis to go 
through. 

10 Short-scale fluctuations 

We now prove Theorem [2 for i < i.e., for Class 3 pairs. These estimates 

will follow from the bounds on NN fluctuations established in Section [6] 
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Lemma 8 There is a constant (3o such that for (3 > (3o , \x — y\ = i < (3^/'^ and 
3m < (3^^^ , we have 

{Bl^) < 2 . (10.1) 

More generally let (xi, j/i), . . . , be n pairs with \xj — yj\ <(forallj, and 

let the interiors of the corresponding C-admissible regions Rx^yi, Rx2y2^ ■ ■ • i ^x„y„ 
be disjoint. Then 

(n^'"^.)^2". (10.2) 

Proof. As in Def. [31 let x be the characteristic function of the interval [0, 1] 
and let (with a parameter S to be defined shortly) 

Xpg^x{{l + 6)-'Bp,) , \p-q\ = l, (10.3) 

and Xpg = l- Xpq ■ Using Xpq < 1 we have 1 < n(pg) Xpq + E(pg) Xpq and 

where the product and the sum are over all nearest neighbor pairs (pq) in R^y ■ 
We estimate the first term on the right-hand side of (|10.4p by applying the 
strategy of the proof of Lemma [5] to show that 

{Bl^ (1 - 3m G,y) n^^^^ Xpq) < 1 ■ (10.5) 
To bound Gxy , note that on the support of Xpq we have 

< ^{tp - tqf < cosh{tp -tq)-l<S . (10.6) 

Thus \tp — tq\< ^/2S and \tz — tj \ < £^/25 for z — x,y and all j G Rxy ■ Now let 
us require 

£^/26 = 1 , or 6^ , (10-7) 

smce I < Thus we have a uniform lower bound on the conductance (17.51) . 

It then follows that < Gxy < G / (3 with G independent of /?, and (jlO.Sp gives 

(^'r IIjp,) < (1 - 3m G/PT' < 3/2 (10.8) 

for f3 large. 

For the second summand of (|10.4p we use 



P/2 

A-pq 



x{{l + 6)-'Bp,)<{^^^ . (10.9) 
The factor -B|™ is estimated by repeated application of ([5.14^ : 

2Bxy<\{ 2Bp^q^ , (10.10) 
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where the product ranges over a set of NN pairs connecting x and y. By com- 
bining (|10.9p and (|10.10[) and then using the result (|6.9p for NN pairs we have 

93m(^-l) , 

/ jD3m \ ^ _f / n/3/2 TT n3m 

\^xy Apq/ - \ PQ LLj Pjl] 



< 



(i-3TO//3) \l~Sm/py 



{1 + S)f>/^ V2 

< e^^VWa (10^^^) 

Since 3m^ < P^/^P^/^ by hypothesis, and (5 = by pTTf)) . we see that the 

expression (|10.1ip is less than exp(/3'^/^ — /?^/^/6). 

Combining our estimates on the two terms on the r.h.s. of (|10.4p we have 

3/2 + 3^V^'^-''^'^/«<2 
for large enough /3. The factor 3£^ < 3/3^/4 

comes from the sum over all NN 

pairs in Ra:y . 

When several disjoint regions are present, the bounds over disjoint regions 
factor, and we can get the same result using the same argument. ■ 



11 Induction hypothesis and some prehminary 
estimates 

The argument in the last section cannot be repeated for all values of £. In order 
to control all scales we need an inductive argument. 

Induction Hypothesis: Let Xi,yi {i — 1, . . . , rii), pj, qj {j — 1, . . . , 712), and 
Tk, Sk (k — 1, . . . , 713) be pairs of Class 1, 2, resp. 3, in the sense of Theorem 3. 
Then the bounds jO]), hold when \x^ - yi\ < i for alii < ni . m 

The induction is on ^ = max^ jx^ — yi|. The said bounds were already estab- 
lished for £ 1 (NN case. Section Ej) and £ < /J^/'* (Section [TO]). Assuming that 
the Induction Hypothesis holds up to scale £, we shall prove (in Section [T2)) that 
it holds up to scale £+1. This will complete the proof of Theorem [4] and, as an 
immediate consequence, Theorem [T] 

The idea of the proof is the same as in Section [TOl If the pair xy is protected 
by a Xxy factor (Class 2), then we apply Lemma [5] in Section [5] 

To get the unconditional estimates we must study the situation when Xxy is 
violated. This violation may happen at any scale from 1 up to £. To quantify 
this we introduce the following definition. 

Definition 5 A point x d A is called n-good if 

Bxy<a\x-y\'' (n.l) 
for all y € A with distance 1 < |a; — < 4" from x. 
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Figure 4: a) In a 3D cube of side 4" we select 8 cubes of side 4'^~^. b) Here we 
see an example of a rooted tree (on a 2D square) with coordination number 5 or 
1 at each vertex, and the corresponding set of subsquares. The root corresponds 
to the large square. 

Definition 6 For a cube Rn of side 4" we define Xr^ to be the indicator func- 
tion of the event that there exists no n-good point in Rn ■ 

Our goal in the present section is to bound the expectation of the indicator 
function Xr„- brief we will achieve this by estimating Xr^ by a sum of 
products of factors of B^y and then using (|9.5p . The details are as follows. 

A 3D cube i?„ of side 4" can be expressed as a union of 4^ disjoint subcubes 
of side 4"~^ . It is clear by inspection of Fig. 2] that we can select 2'^ = 8 of these 
subcubes, say Ri,_i {i = 1,...,8), so that dist{Rl_i, R{_i) > 4""^ {i ^ j). 
Our approach now rests on the following simple observation: if there is no n- 
good point in i?„ , then there is either no (n — l)-good point in any of the 8 
subcubes R^^i , or else there exists at least one bad pair {x,y) e i?„ x A at 
scale 4"~^ < \x — y\ < 4". Thus, Xr^ is bounded by the inequality 

8 

xk<5L+nx<_,> (11-2) 

2=1 

where 

sk.= E ^s,- (11-3) 

X e Rn, y e A 
4"-i < |a; - y| < 4" 

We iterate (|11.2p inside each cube -R^„i , thus selecting 8^ subcubes of side 
4"^^, and we keep repeating this procedure until we reach cubes of side 4° = 1 
(i.e. points). We denote by Rn-k [k = 0, . . . , n) the set of 8*^ cubes of side 
4""'^ obtained in this way. i?„ = is the starting cube. Moreover let 

R — ^k^gRn-k- In this way we bound Xr by a positive sum of products of 
Xxy 1 which in turn arc bounded by Bjp / {a |a; — yj")™. The resulting expression 
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can be organized as a sum over rooted trees picturing the hierarchy of inclusion 
relations of the subcubes. 

The following set of definitions serves to prepare the statement of Lemma [H] 
below. Let V be an abstract set of vertices such that \V\ — \R\. We associate 
by a fixed bijective map each vertex w in y to a cube i?„ € R. We denote by ky 
the scale of the corresponding cube: Ry G Rn-k^- The vertex associated with 
the largest cube Rn is denoted by r (root). Finally let A{v) (ancestor of v) be 
the unique vertex in V such that Ry C Ra{v) (see vertex 1 and 1' in Fig. \^). 

With these definitions we can introduce 7^ the set of labelled rooted trees 
on some subset of V with root r, such that the root has coordination number 
dr = 8 OT dr = (in which case the tree is reduced to a single vertex) and 
the other vertices have coordination number dy — 9 or dy ~ 1. Moreover if v 
belongs to the tree then there must be a tree line connecting v to its ancestor 
A{v). The maximal distance of a vertex from the root is n. Let Lt denote the 
set of vertices in T with dy — 1 (the leaves) or rf,- = (then Lt contains only 
the root). Let 14 be the set of vertices in T at distance k from the root. With 
these definitions the tree is completely fixed by the leaves Lt (or equivalently 
by the choice of the coordination numbers for each vertex). See Fig. |4)d for an 
example in the case of dy — 5 instead of 9. 

Lemma 9 With the definitions above we have the inequality 

Ter„ veLT 

where Ry is a cube of side 4"^''" and Sfj^ is defined as in (|11.3p . with n replaced 
by n — dy . 

Proof. Our trees T e 7^ are constructed by iterating (|11.2p . In each iteration 
we get to choose between the first and second term of the r.h.s. of (|11.2[) . 

The construction starts with the root of the tree. In the first step of the 
iterative scheme, if we pick the first term of ()11.2|) then the construction ends 
and we have produced nothing but the trivial tree (the root). If we pick the 
second term, we have a product of 8 different indicator functions Xr ? one for 
each subcube in Rn-i ■ We represent them in the tree by attaching 8 vertices 
to the root. Each vertex is then associated to a subcube by lexicographical 
order (see Fig. IDd). Now we repeat this procedure at the end of each branch 
and, continuing in this way, construct a tree. Whenever we pick the first term 
in (|11.2p . the corresponding branch of the tree terminates and we produce a 
terminal vertex or leaf. If we pick the second term, we generate a vertex of 
coordination number 9. The iteration stops when we reach scale n. m 

We are now in a position to state and prove the main result of this section. 

Proposition 1 Let d = 'i and let the parameters a, to, a be chosen such that 
ma > 4:d and a > 2'^"+('*+^)/™. Assume that the induction hypothesis (|9.5p 
holds up to scale i = 4"+^. Then 

{xrJ < 2-("+i)"™ , (11.5) 
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and if Rn(k) , k — 1, . . . , N, denotes a family of cubes of side 4"('^-' < £ such that 
dist(i?„(fe),i?„(fc,)) > max(4"W,4"('=')) (11.6) 

then 

Proof. Note that the average of (|11.3p is bounded by 

(^«.> < E 

X e Rn, y e A 
4"-i < |a;-y| < 4" 
2 

for all n > 1. For n = the sum is over e A at distance |x — y| = 1 and we 
obtain the same bound (see below). 

First, consider the minimal scale n = . In this case i?„ = Ra is just a single 
point, the sum over y is a sum over 2d nearest neighbors, and we simply have 

{x'kJ<{SnJ<2d^<2—. 

Thus (fTT31l holds for n = 0. Now let n > 1. By (flT2l) we have 



Assuming that (x^ _i) — 2~""™ holds, we conclude that 

Hence the result is true for all n by induction. 

To prove (|11.7p we apply (|11.4p to each term xf? • Then we can again apply 
(|9.5p as long as the corresponding diamonds are disjoint - this is ensured by 
the procedure for choosing subcubes and by the constraint dist(i?„(fc) , Rn(k')) > 
max(4"*^''\ 4"^'^ )). This concludes the proof of Proposition [TJ ■ 

12 Proof of the Induction Hypothesis 

In this section we shall establish the induction hypothesis of Section [11] at scale 
£ assuming that it holds up to scale £ — 1. Since our regions RxiVi are disjoint 
by assumption, we will be able to re-express each factor Bx^y^ as a sum over 
non-overlapping regions where our induction hypothesis applies. To simplify 
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the notation let Xi = x and yi — y. We will assume that £ — \x — y\ is large, 
i.e., 1^1 > (The case of small £ < (3^/* was dealt with in Section[lOl) 

For z = x,y we recall the meaning of the regions R^y from Def. [21 Section[71 
To ensure that the new regions produced by the analysis below remain inside 
the original region R^y we need to introduce the following subsets. 

Definition 7 Let R^y he a diamond region as described in Def. Then we 
define R^y for z = x,y as 

Rly = {j e Rly : ^{jz,xy) < 7:/8 for \j - z\ > 10} , (12.1) 

where Z(jz,xy) is the angle between the lines jz andxy. 

This definition roughly selects (at distances larger than 10) a double cone which 
is obtained by rotating around xy a 2D diamond with vertices on x and y and 
opening angle 9 = tt/8 (see Fig. [5^). The condition \j — fc| > 10 ensures that 
R^y U {z} is connected. 
We also define 

= n n (12-2) 

for 

Xxj J Xyj defined in (18. 1|) . Note that R%y U R^y U {x, y} is a 5-admissible 
region in the sense of Def. [T] (Section [7]), so Lemma [5] and Lemma [6] can be 
applied to give 

{B^:^yU.,y)<{l-mC/pyK (12.3) 

Now let x'xj = 1 ^ Xxj and Xyj = 1 ~ Xyj ■ The next lemma is nothing but 
a combinatorial identity based on the following partitions of unity: 

1 - n (x-j- + ^Sj) ' 1 = n ^^y^ + ' (12.4) 

j<^R%y j<^Rlv 

where R%y and are the regions defined above. 
Lemma 10 The identity function can be written as 

i=u.j,y+ ^ xlb n + n ^2/'' n n ^y^ ■ 

b£R%y j,\j-x\<\b-x\ b^Rly i&R^y j,\j~y\<\l'-v\ 

(12.5) 

Proof. We start from (|12.4p and expand the first product over j, beginning 
with small \j — x\. For each factor Xxj + Xxj have two possibilities: either 
we pick Xxj 7 in which case we proceed to the next factor and repeat, or else we 
pick x'xj and then we stop expanding and leave all the other factors (with larger 
|j — x|) in summed form Xxj + Xxj = 1. In the resulting sum there is the term 
n,pRx Xxj ■ This we multiply by the other product (over j e Rly) in (|12.4p . 
which we expand in the same way. 
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In total, we have either picked a factor % for all j G R^y and j G i?^^ (this 
results in the term Uxy), or we have picked a term Xxj or Xyj somewhere during 
the course of the expansion process (this gives all the other terms). The point 
where we stopped is denoted by b (where b stands for 'bad') because Xxb > or 
Xyt > means that there is a large deviation at that point. ■ 

Using the equality (|12.5p of Lemma [TU] we can rewrite BJ^y as 

=i?;>,, + 7^(x,y), (12.6) 

where 7^(2:, y) is defined as 

beRly 

We have Uxy) < (1 — toC//3)~^ by (|12.3|) . without any need for an inductive 
argument. Thus if we can prove that {TZ{x,y)) < 2 — (1 — mC / f3)~^ our proof 
will be complete. The desired statement is formulated in the next lemma. 

Lemma 11 For large (3 the remainder ()12.7p is bounded in average by 

{n(x,y)) < 1/2. (12.8) 

Moreover, TZ(x,y) can be written as a sum over products ofBx'y' with \x' — y'\ < 
i — 1 in such a way that the corresponding regions Rx'y' are disjoint. 

Now, using this lemma and arranging for (1 — mC/P)'^ not to exceed 3/2, 
we have 

(B™ ) < 3/2 + 1/2 = 2, (12.9) 

thus completing the proof of the Induction Hypothesis for the case ni = 1 and 
n2 , ri3 = . The general case - (|9.5p - is done in the same way. 

Before starting the proof of Lemma [Tl] we provide some orientation and 
motivation. To bound the expectation of Xxb '^^ could try to use (|5.14p : 

< 2"^B2Bry , (12.10) 

while Xxb — 1 ~ Xxb can be bounded by 

x'xb<B'^,a-P\b-x\-"P (12.11) 

where we used Def. [3] of Section [H] Since both |& — a:| and |fc — y| are smaller 
than i it is natural to try to apply the Induction Hypothesis. However we face 
at least two problems: 
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[A] . The Induction Hypothesis does not cover when p > . Indeed, a 
factor Xxb — n Xxj Xbj rnust be present in order for (|9.4[) to apply in such a 
case with < p < 2m. Notice, however, that while we have no immediate control 
without the missing factors Xbj, the factors Xxj for j € R^y and |j — < |& — a;| 
are already in place. To overcome the problem, we shall introduce the needed 
factors Y[ Xbj by the same partition of unity scheme that was used above. 

Before embarking on that scheme, let us quickly evaluate the situation which 
emerges after insertion of Xxb ■ We can then choose p = 2m , the induction (|9.4p 
applies, and we get a small contribution 

{B^yXlbXxb) < 2"^a-'P\b-x\-^^{B":+''XxbBry) 

<4.2ma-2m|^_^|-2am (12.12) 

for \b — x\ large. Here we used (j9.4l) since xb is of Class 2 and by is of Class 1. 
Note that the expression (|12.12p is summable in b for ma large. Moreover, the 
factor a"^™ ensures that also the contributions for l^* — 2^1 = 0{1) are small. 

[B] . The second problem is that, in order for (|9.4p to apply we must make sure 
that we can find inside Rxy two non-overlapping regions R^b and Ri,y of which 
the former is C-admissible and the latter of diamond type. Moreover, since we 
have Y[ Xxj only for j G R^y we must ensure that R^i, is inside the reduced 

region R^y U R}jly U {x, y} in R^y ■ Since this might have to be repeated many 
times at smaller and smaller scales, we must be sure that all regions remain 
(5-admissible (or, put differently, we do not want 6 to be scale-dependent). We 
will see in the next lemma that this can be arranged. 

Lemma 12 For a diamond Rxy consider the subsets R%y and R!^y of ()7.4|) . 

[1] Let w be any point in R^y . Then we can always find a point a € Rxy 
and regions Rxw , Rwa , and Ray , such that Rwa and Ray are diamonds inside 
Rxy , Rxw inside R^y U {x} and is C-admissible, and the three regions have 
disjoint interiors (see Fig.\^. 

[2] Let w be any point in R^y such that \w — x\ > (xw not of Class 3). 
Let W1W2 be a pair in R^y U R^y U {x, y} such that \wi — w\ < \w — and 
\w — x\^^'^ < — W2I < \w — x\/5 (see (112. 18p ). Let R^^^^^ be the corresponding 
diamond region. Then we can always find ^ points g Rxy (* = 1, • ■ • , 4j, such 
that all of the regions Rxai, Rajaj+i (j ~ l,---,3j, and Ra^y , are diamonds 
with disjoint interiors and do not overlap with RwiW2 (see Fig.\d^. The same 
can be done for x € R%y . 

Proof. The most dangerous situations are shown in Figures [5] and [51 It is a 
simple geometrical argument to see that the region Rxw in Fig.[5]is C-admissible, 
as the angles 9i and 82 are never smaller than 7r/8 , see (|12.ip . In the cases shown 
in Figs. [5^ and [5)3 one has to check that the diamonds do not transgress Rxy ■ 
This never happens since \wi — a;| < |x — y\/\/2 and \wi ~ W2\ < |wi — a;|/5. ■ 
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Figure 5: Wc need one intermediate point a. The two angles Oi and 62 are never 
smaller than tt/S. 
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Proof of Lemma llll We split the sum over bad points b in (|12.7p into several 
groups of terms. 

Case 1. The bad point b is located close to x, i.e., |6 — < /?^/^. Then we 
can bound by (|5.14p and x%b (|12.1ip . which gives 

BTyXlb < ^-^a-Plb-xl-'^PEZ'^Br^. (12.13) 

To apply the Induction Hypothesis we need to select inside the diamond Rxy 
two regions R^t and Rty ■ The first one, Rxb , need only be C-admissible (since 
xb is of Class 3, see Section [3]), so it may be a deformed diamond (Fig. [5)3). 
On the other hand, Rby has to be diamond-shaped, since \b~y\> (Class 
1). To make the requirement of diamond shape conform with our constrained 
geometry, we must add an intermediate point a as in Fig. [5] with w = b, and use 

Bll < r'BZB:^^ . (12.14) 

We have seen in Lemma [T^] that we can always find such a point a, so the 
induction (|9.4[) does apply. Note that since |6 — a;| < (3^^^ there will be no 
additional induction on R^b ■ Therefore there is no risk that the region might 
get more and more deformed by the induction steps and (5-admissibility might 
finally be lost. Thus we have 

B^yXlb < 2'"'a-P\b - xl-^^BZ-^'BraBTy ■ (12.15) 

The situation for b near y is analogous. Summing the contributions from b near 
a; or y we obtain 

z=x,y |6-2|</3i/4 ^ ' 
< ^ ^ 4"a-2m|5_^|-2am22(i^_^) 

z=x,y |b-z|</3i/4 

4™ n 

<2\l + p)—K^ l^-^l'~'""<(4/a'r^;<^, (12.16) 

|b-a:|=l 

where in the second line we used ()9.4p and p = 2m . We can accommodate 
m + p = 3m > m without any protection factor x since bx is Class 3. In the 
third line, Ki\b — x\'^ is the entropy factor for the 3D sum over bad points at 
distance \b — x\, the factor K[ is a constant of order unity, and we used that 
4/a^ < 1 and m > 4d/a is large. We bounded the expression by p/10 for 
convenience; since both a and m are large, the factor {A/a^)™K[ is in fact very 
small. 
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Case 2. The first bad point b is far from x (i.e., |6 — > (3^^^) and also far 
from y. Let us consider the case b € R^y for definiteness. (The other case, 
6 G R^y , is treated in the same way.) Again, we have to estimate 

(sTyXlbU ,. ,X.,\ • (12.17) 

As was observed above, if we succeeded in promoting the last product in the 
average to a complete factor 7 then we could apply the Induction Hypothesis 
as in p2.12p . In order to satisfy the hypothesis of LemmaO Eq. (|7.6p , we should 
have a constraint Xjh for all \i — b\ < \b — x\/^/2. Actually, from the remark 
after the proof of that lemma we only require Xjb for \j — b\ < \b — x\/5 since 
we know that all — tj fluctuations are good up to |j — ~ |6 — 

Guided by the idea of partition of unity (cf. (|12.4|) - (|12.5p ). we will first check 
whether there is some large fluctuation > at large scale near b. If no such 
event occurs, we proceed to the step of checking fluctuations at intermediate 
distance scales. Then either all intermediate distance fluctuations are good 
too (and we have the desired factor x), or there must be some bad event at 
intermediate scale. In this last case we will see that many bad events must 
happen. We will now make this more precise. 

Case 2a. The nearest bad point b is far from x (and y), |6 — a;| > P^^^, and 
there is a large scale bad event near b. This means that Bj/^ > a |j — for 
some pair j, k e R^y U R'^y such that 

\j - b\ < \b - xl^/"^ and \b-x\^/'^ <\j -k\<\b-x\/5 . (12.18) 

Now, using (|5?T4l) and (|12.1ip . 

3 

BTyX-jk < 2''"i?ra, n Ka.^BZyB^Jl \j - a"™ . (12.19) 

To apply the Induction Hypothesis the corresponding regions must all be 
diamonds (all pairs are Class 1). By the assumptions made on the pair jk, 
Lemma [T^ guarantees that we can choose the four intermediate points S Rxy 
(i = 1, . . . , 4) so that aU of the regions i?a;ai , Ra^a^+i (j = 1, ■ • • , 3), and Ra^y , 
are diamonds with disjoint interiors and do not overlap with Rjk (see Fig. [5]). 
Since the regions are non-overlapping and £ > \j — k\ > |6 — our induction 

hypothesis yields 

{B™X,\) < 24™26|6-a;r"'"/2a-™ . 

For large m the value of the sum over b is small. 

To estimate the entropy factor, note that there are less than |5— a;|'^+'^/^ pairs 
jk satisfying (|12.18p . Altogether then, the present partial sum of contributions 
from r = \b ~ x\ > [3^^^ is bounded by 

(4/a)™7^2 r('^-i)+'^+'^/2-"'"/2 ^ o{(3-^/^) < — ■ 
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Note that there is nothing special or optimal about the exponent 1/4 of 1//3 - 
it is just convenient. 



Case 2b. We now suppose that |& — a;| > (3'^^^ and there is no large deviation 
near h, i.e., Bjk < a |j — fc|" holds for all j, k subject to (|12.18[) . This implies 
that at long scales \j — k\ > \b ~ xj^/^ we have Xjk — 1- It remains to check 
whether Xjfe holds also at shorter scales |j — fc|<|6 — a;|^/^. 

First we consider the case of there being a point g {g stands for good) in R^y 
with \g — h\ < \b—x\^/^ such that Xgh = 1 holds for all h with \g — h\ < \b~x\^/^. 
We then have in particular that Xgb — 1, and so by Def. [3] 

Bgt < a |& - g|" < a {\b - x]^/^)" ^a\b- x\°'l'^ . 

This inequality combined with the constraint — 1 and (|5.14p yields 

Bxh a\b — xY^ ,„/9 
Bgh a\b-x\°'!^ 

Thus we have 

xSfcX.s < '^-"\B^^x.g)B'^yX'^^ < 2^"^{Bl^x.g)B^y \b - x\-^^ . (12.20) 

Now we have to be somewhat careful about the choice of the regions Rxg and 
Rgy , as they may not have the canonical diamond shape. For Rxg this is not 
a problem, because of the presence of Xxg {xg is of Class 2). All we need is 
that Rxg be C-admissible. On the other hand, Rgy is (as in Case 1) slightly 
more delicate. To be sure that we deal with diamond-shaped regions, we add 
an intermediate point a as in Fig. [5] and use B™ < 2™i?™i?™ . We have seen 
in Lemma [T^ that it is always possible to find such a point a. 

The regions Rga and Ray are of diamond type, so induction applies. Note 
that Rxg comes with a Xxg factor {xg is of Class 2) so no additional induction 
is required for it. Therefore, as in Case 1, there is no risk that the region might 
get more and more deformed by the induction steps. It should be emphasized, 
however, that Xxg is not exactly the same as in (jS.ip . but rather is given by 

Xxg = n n • 

-RxH^J ■ \j-x\<\b-x\ R^yBj ■■ \j-a\<\b-x\/5 

Since — 5I < \b — x\^^^ and |& — x| > (3^^^ we have |(7 — ~ \b — x\ up to 
a correction factor of order 0{\b — x\^^/^) < 0(/3~^/®) <C 1. Therefore Xxg is 
equivalent to the following constraints: 

Vj e Rxg , \j -x\< \g - x\fi : Bxj < a\j - x]" , 
and yj e Rxg : \j - g\ <\g - x\f2 : Bg^ < a\j - g\" , 

with /2 = 1/5 and /i = 1 — 0{l3^^^^). From Remark 7.3 we know that Lemma 
[S]and hence Lemma [5] still hold, so we can apply the induction and 

(BTyX'xbXxg) < 2^™(1 + p) 2' \b — . 
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There are 0{\b — a;|'^/^) choices for g, so the sum over these contributions is 
bounded by 

2^"'K3 J2 r''-i+'^/2-am ^ 0(^-1/4) < A . (12.21) 

r>/3i/4 

Case 2c. The last case to consider is the situation where no such point g 
exists. In that case we can always find a cube i?„ C R^y U K^y which contains 
the point b and has side 4" = \b — x\^^^ such that i?„ is at least at distance 
^1/4 < 4" ^ £ from the boundary of Rxy and contains no n-good point (see 
Def.Oin Section [TT])- Then by pTI]) and (|12.1ip with p = m, we have 



Xr, 



< 



TT , c ^ V- TT 



E n.„*usE i: n 



where according to (|11.3p the sum over configurations of pairs {jy, ky) € Ry x A 
is constrained by 4"""^ < \jy ~ ky\ < 4"", with Uy the scale of the leaf v. 
Since all cubes Ry are inside the small region i?„ and all pairs jy , ky satisfy the 
conditions of Lemma 1 121 for the pair wiW2, we can proceed as in Case 2b and 
select 4 intermediate points aj, j = 1, . . . ,4 such that the corresponding regions 
are diamonds and do not overlap with any Rj^,k^ (see Fig. [S]). Then by (|5.14p 
we have 

B.yXR„<^ B.a,---Ba,y I.™ ' (12.22) 

and we can apply the Induction Hypothesis. By Proposition [T] in Section [TT] we 
have 

(B^XrJ < 24"2-""™ , n^ln\b-x\> \n{(3^^'^) . (12.23) 
Therefore we have enough decay to control the entropy factors: 

r>/3i/4 

This concludes the proof of ^2M- From Eqs. (|12.15p and (|12.19p - (|12.22p . 
we see that TZ{x, y) can be written as a sum over products of such B^'y' with 
|a^' — y'l < i— 1 that the corresponding regions Rx'y' are disjoint. This concludes 
the proof of Lemma fTTl and the Induction Hypothesis. ■ 



13 Proof of Theorem 2 

Now that we have estimated (-B^*^) for all x,y we need to estimate (cosh^(ta;)) 
for moderate values of p < 10. If we suppose that the field t is pinned at some 
point jo , so that tj„ — , then Theorem 2 follows directly from Theorem 1: 

(cosh^'t^) = (coshf (i^ - tj„)) < 2 , 
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for any x in the lattice (since Theorem 1 does not require bounds on e). When 
the field is not pinned, we need £ > and some conditions on the volume. The 
rest of this section is devoted to this case. As in the proof of Theorem[T]we will 
first prove bounds on conditional expectations. 

Definition 8 A point x Cz A is called 'good at all scales ' if 

yjeA\{x} : B,j<a\j^x\" 

(see also Def. in Section\7l]). The corresponding characteristic function is 

where the factors Xxj ore those of Def. (Section\^. 
Lemma 13 Let x be good at all scales, and let 

= cosht^ + ^e*" si . (13.2) 

If /3 ^ 1 and e > Spa L^'*+", then for any < p < 0{/3) we have 

{Blxx)<2. (13.3) 

Proof. The proof uses a combination of ideas already present in the proofs of 
Lemma[S](Section[71) and Lemma|n](Section|S]). By supersymmetry (Proposition 
m Appendix [C)) we have 



LjGA\{x} - 

where Zx is defined in (|2.11|) and Xxj (|8.3p . Following exactly the same steps 
as in the proof of Lemma [6l we obtain the inequality 

l>{Blxx{l-pGx)) , Gx = ^[Sx;D0,,{t)-'Sx], (13.5) 

if pGx < 1- We must now bound the Green's function Gx using the constraint 
Xx (as we did in Lemma[5|). For this purpose define D ~ e^*^^BxDi3,^{t) by 

[v-bv]^e'-BxPy^, e*'+*^-2*^(«,-t;,)2+eS, V, e*''-*^^^ , 

^ ^ [ij] ^ — 

and note that Bxj < a\j — implies the bound 

e*^-*- > {2a\j-xn-^ . 
We then follow the proof of Lemma [5] and introduce a telescopic sum 

N 

Sx = {Sx - h) + (/i - /2) + . . . + {In-1 -In)+In = J2 ' (13-6) 

n=0 
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where J„ is the (normahzed) indicator function of a cube of center x and side 
2", and pn — In — In+i ■ There is no need to introduce / as we did in the 
proof Lemma [51 as we are now working not on R^y but on the whole volume. 
The sum terminates on reaching the system size 2^ . Note that for n < A'' we 
have Y.jPn{i) = and \\pn\\l < 2""'^ = 2-3". The function p„ for n iV is 

constant: pnU) — InU) — \M^^ ^'^^ J ^ ^■ 
Now, by the Cauchy-Schwarz inequality, 

[<^.;^-'<5.] < fe[p„;^-Vn]'/') . (13.7) 

For n < iV we use the bound [pn;D~^Pn] < II (-D""^)p„ || , where || (I)"^)p„ || < 

2^"+^""ci . Then for the sum of terms with n < N we have 

N-l N-l 

EK;^-V.]^/^ < ^ E V2-"<^-^-^"^ < ^ (13.8) 

uniformly in N since d = 3 and 2q; ^ 1. For n = A'^, on the other hand, we no 
longer have orthogonality to the constant functions ('zero mode') and therefore 
must take recourse to the e-term in D : 

[pn;D-'pn] < (e|A|mine*^-*^)"' < 2a-^ < ^ . (13.9) 

Hence < (7i/\//3 + VV^)^ < l/(2p) for /3 > 1. So, 

|(i?Sx.> < (S?X.(1-PG,)) <1 (13.10) 

by (|13.5p . and the lemma is proved. ■ 
With this lemma we can finally complete the proof of Theorem [2l i.e. the 
bound on the unconditional expectation of cosh^ . 

Proof of Theorem [2l We recall from Dcf. [5] (Section [TTjl that a point x is 
said to be n-good if B^j < a\j — x]" for all j G A subject to 1 < \j — a;| < 4". 
A point X is good at all scales if B^j < a\j — x|" for all j e A \ {x}; we then 
say that x is A^-good. 

We proceed as in Lemma [TT] (Section [T^ : 

(cosh'' = (xx cosh" + {xl cosh" t,) , (13.11) 

where Xx ensures that the point x is iV-good. Then by Lemma fT3l we have 

(Xx coshP < 2 . (13.12) 

It remains to estimate the second term, {x^coslf tx) ■ We prove in Lemma fT4l 
below that this term is bounded by a constant. Once this has been accomplished, 
the proof of Theorem 2 will be finished. ■ 

Lemma 14 Let Xx — ^ ~ Xx , with Xx defined by ()13.1|) . Let (3^1 and 
e> 8 • 4 • 10 a L°'-'^. Then for any < p < 10 we have 

ixlcoshPtx) <l/2 . (13.13) 
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Proof. While x% means that x is not iV-good, it is still possible for other 
points in A to be A''-good. If g ^ a; is the nearest such point (as seen from x), 
then none of the points inside the ball K?_. of radius \g — x\ and center x is 
A^-good. Denoting the indicator function for the latter event by x'k'^ have 
the identity 

E ^9XK^^_^+X{xh (13.14) 

where the last term accounts for the possibility that there is no A^-good point 
in A at all. Thus we obtain the decomposition 

ixl coshP t,) - {^s Xk^^_, cosh^' t,) + (J{^ X ■ cosh^' t^) . (13.15) 
We will prove that both of these two terms are bounded by 1/4. 



1. We consider the first sum. Using coshtj; < 2cosh(ta; — tg)coshtg and ap- 
plying the Cauchy-Schwarz inequality twice, we obtain 

{xK^^_^^XgCosii''tcc) < (x^x_^^Xs)^^^XsCOsh^^'t^)^/^ 
^ (^^-.,^3 )'^' 2^'(x, cosh2^'(t, - tg) cosh^'' Q'^' 
< 2^ (X^._,X.>'/'(X. cosh^^ cosh^^(^. - ^.))'/' 

where in the last step we used Lemma [13] and Theorem [1] and we introduced 
Co = 1 and Cp = 2^/^ for p > I. It remains to bound 

g n>0 4"<|g-a:|<4" + i 

Let be the cube centered at x of side 4". Now fixing a point g with 4" < 
Is' — a;| < 4"+^ we have i?fj C K^g_^^ , and we distinguish between two cases: 

la. The interior of K^g_^^ is void not only of A^-good points but also of n-good 
points. Let Xn denote the corresponding indicator function. Then for n > 1, 
using Proposition [1] we have 

(X^r ,X.X^)<(x^s)<2-("+^)"", (13.17) 

^ \g — x\ ' ^ Ti ' 

where Xr^ is given in Def . [51 Section [TlJ For n — the cube Rq contains only 
the point x, so 



|2 — a:| — 1 — a;] — 1 
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lb. There is at least one n-good point y inside . Let Xn{y) be the 

corresponding indicator function. Because the point y cannot be iV-good, there 
must be a first scale q > n so that y is g-bad. Thus there exists a first point b 
at distance | & — 2/| > 4" with Byi, > a | & — y|". It follows that 

E {XaXK^^_^^Xniy)) < J2 E a"M6-V"' 

And I. 9t-'^^11-„ 

r>4" 

where the factor 4"'^fci comes from the sum over y and r'^^^k2 comes from the 
sum over b. Inserting these results into ()13.16p we obtain 

ri>0 4"<|g-a:|<4'»+i " n>0 

where 4*^"+^)'*fc3 comes from the sum over g. This will be no greater than 1/4 
provided that am is large enough. 

2. To complete the proof, we have to estimate the last term (Ylj x'j (^osh^ t^) 
in (|13.15p . By Proposition [T] the probability for no A^-good point to be found 
in a cube A of side L = 4^ is bounded by 2'^"™ = i-am/2_ n^^^e 

(n^.X^osh^t.) <{cosh'n,y/\Y[^x'jy^' < (cosh2Pi.)^/'L— /4 . 

To get a bound on the expected value of cosh^'' we once again use supersym- 
metry (Proposition [21) , as follows: 

e'^- = (e^"^) = (e'^-^^(l-7£G^)), 

where we choose < 7 < 1/2, and G'^ — e*=^' [Sx ; I?/3,e(i)~^(5a:]. Since the operator 
Df3,e{t) — se^^^xiSx ; •] is non-negative, by Lemma[T]we have eG^ < 1, so 

(e7^(B.-i)) < (l-7)-i . 
Also, cosh^^'ia; < (2p)! (72)"^^ e''''^^=" by an elementary computation, and hence 

(cosh^Pt,) < (2p)! (7£)^2p^g7es.^ < o{e-^P) . 
We thus finally obtain 

iYl.X'jCOshPtx) < (cOsh2Pt^)l/2i-am/4 < 0(^-p)^-am/4 ^ -^^4 ^ (13.I8) 

since am is large and e > L"^'^. This concludes the proof of Lemma fT4l ■ 
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Remark. In the proof of Theorem [5] the e term (zero mode) appears only in 
two places: (|13.9p of LemmafTSl fthe last term in the telescopic sum) and (|13.18p 
(when no TV good point is present). The inequality (|13.9p is the reason why we 
cannot take e = 0{L^'^) but must take e — 0{L°'^'^). 

14 Proof of Theorem [3] 

Finally we can prove the bound on the Green's function Cxy of ()1.9|1 . Let / be 
such that /(j) > for all j & h-. We need to estimate 

[/; Cf] = ([eV; Df^^Sr^e'f]) = {[W; GtW]) , (14.1) 

where Dp^e{t)-^ = Gt was defined in dTT]), and W{i) = e*^/(j). 

14.1 Upper bound 

Let Lq = —(31s. + £ and Go = L^"^ (as defined in the statement of the Theorem). 
Now 

[W; GtW] = [LfiGaW; GtW] = (3 [V{G^W)-V{GtW)] + e [GqW; GtW] 

= PJ2 ^n'{GoW)V,,,{GtW) +eJ2iGoW)U){GtW)U) 

in') i 
„ a^ ^n'iGoW) VjriGtW) f {GoW)m f iGtW){j) \ 

Since \a ■ b + c ■ d\ < {a ■ a + c ■ cy/^{b -b + d- df/"^ we have 

V to') i I 

TherefortQ 

in') j 

Now 

|V,y(GoVK)| <^|(Go(j,fc)-Go(/,A:))|iy(fc) <const^i/,,M/(fc) , (14.3) 

k k 

where we defined H^k = /3~^(|j-/cp + l)~ie-^~IJ-'=l, £ = (e/2/3)i/^ and we used 
|(Go(j- fc)-Go(j',fc))| < const i/jfc . 



^We thank S.R.S. Varadhan for explaining the inequality 1114. 2| l to us. 
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By inserting into we get 

[/;C/]< const /3 ^ /(fc)/(0 (e(*'=+*-*^-*.')) 

+ eJ2GoU,k)Go{j,l)f{k)f{l){e'^+'^-'^) . (14.4) 

By Theorems [T] and [5] the expectation over the field t is uniformly bounded. 
Now we can sum over j : 

^ HjkHji < const Go(fc, 0, J2 ^"(j' ^) = 0, (14.5) 

j 3 

where Go = (-/3A + e/2)"^ Note that Go < Gq . We finally obtain 

[/; Cf] < const [/; Go/] + e[f; G^/] < 2 const [/; Go/] • (14.6) 
This completes our proof of the upper bound. ■ 

14.2 Lower bound 

Let Xx be the characteristic function ensuring that a: G A is good at all scales 
(see (|13.ip ). Recall that if > then 

e*^-*=" > (2a(l + |j-2:|"))-i (14.7) 

for all j € A. We have the inequality 1 = Xa; + — Xa; • Inserting it into (|14.ip 
we obtain 

([M^;Gtiy]) >(x.[W^;GtM^]) -^(x.W^(j-)W^(fc)Gt(j-,A;)) (14.8) 

jk 

where f{j) = (1 + |j - x\-)-^f{j) and D'^ = e^'^Gt = {e-^*- D0,,{t))-\ In 
the first line we used the fact that Gt is positive as a quadratic form for each 
configuration of t. In the second line we used the fact that this is a sum of 
positive terms since W{j) > and Gt is pointwise positive. Furthermore, we 
applied (|14.7p to estimate W{j). Now, 

(x.[/; A^Vl) = (x.) E([/; A" Vl) > (x.) [/; E(A)- VI (14.9) 

where 

E(.) - ^^^'^ 



iXx) 

is a probability distribution and we used Jensen's inequality. 
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In order to complete the proof we need to estimate (xx) and (^XxDt). From 
Lemma [14] (in the previous section) with p = we know that 

iXx) = 1 - iXx) > 1/2 . 
Moreover {x^D^ < (5^ as a quadratic form and for any function u we have 

[«; {D) u] = - u{k)f (e*^+*-2*^) + eY,u{jf {e'^'^'^) 

Uk) j 

</3c^(7.(j)-"(fc))'+£c'^u(i)2 <ci[w;Gou] , (14.10) 

where we apphed Theorems [1] and El and ci = sup{c, c'}. Thus (D) < ciGo ■ 
By applying these relations we see that E(Z)f) < 2ciGq and hence 

[/;C/]>^[/;Go/]. (14.11) 
This concludes the proof of Theorem [31 ■ 

Remark. If W did not depend on t we would have the quadratic form estimate 
[W; GoW] ci < {[W- GtW]) < C2 [W- G^W] 

with 

C2 = sup ((e-*^-*^') , (e-*-)) , C2 = sup ((e*^+*^') , (e*'=)) . 

Uj),k {jj'),k 

The upper bound follows directly from (|14.2|) . the lower bound from Jensen's 
inequality. 



Appendices 
A Minimum of the effective action 

Let J i~> = f e M (for all j G A) be a constant field configuration. Evaluating 
the statistical weight function on it we get 

g-£|A|(cosht-l) j3g^l/2(_^^_^gg-t) (^ -^^ 

Let t* be the number that maximizes this statistical weight. The condition for 
the first derivative to vanish at t* is 

2sinhr = -e"**Go(a;,a:) (a; g A) , (A.2) 
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where Go = (— /3A + ee * ) > 0. Equivalently, 1 — e^* — Go{x,x), and since 
Go is non-negative, it follows that t* < . We thus infer that 

< 1 -e^*' = Go(a;,a;) < 1 . (A.3) 

Next, we show that the constant field t* maximizes the integrand over the 
full set of all field configurations t — {tj}. For this, we recall the definition (|f .2p 
of the effective action or free energy F^.^ in combination with (|f .lOp : 

= /3V^.. fcOsh(t,-t,)-f)+£V (coshtfc-f) 

^ In Dct^/^ {-I3A + l3V{t) + e c"*) . 
Now we introduce A Gl^^{pV{t) + e (e"* - e"**")) G^^ and write 

Dct(-/3A + /3V^(t) +ec-*) = Det(Go Det(Id + A) . 
Using lnDet(ld + A) < Tr A we then obtain 

Fi3,e{t) > - lnDet(Go ^/^) + e|A| (cosht* - f ) 

+ e^^fe (coshtfe - cosht* - ^(6"**= - e"*') Go(a;, x)) . 

The second line of the r.h.s. is non-negative by ^ J2j ^jj — J2{ij)i^o^^i^i~^j)~^) 
and Go(a;,a;) < 1, and so is the third line by the identity Go{x,x) = 1 — e^* 
and a trivial computation. This proves that Ffj,^{t) is bounded from below by 
Fp^^{t*) = -lnDet(Go^/^) +e|A| (coshi* - 1). 

B Hyperbolic symmetry 

In Section 0] we explained that the H^'^ nonlinear sigma model in the limit of 
vanishing regularization e — > 0-1- acquires a global symmetry by the Lorentz 
group S0(1, 2). We will now exhibit the Ward identities due to this Lorentzian 
symmetry S0(1, 2). (Consequences due to the supersymmetries of model will be 
explored in Appendix lO) To prepare the discussion, the reader is invited to re- 
call the expressions (j2.7p for the functions x, y, ^, rj in horospherical coordinates. 
We also recall that z = coshi -I- e*(|s^ -|- '0V')- 

We now seek the first-order differential operator, Li, generating Lorentz 
boosts in the zx-plane, i.e., 

Li z = x , Lix — z , Liy — Li(^ — Lir/ = . 

It is easy to verify that the unique operator with these properties is 

Li = dt- i^d^ - ipd^ - sds . (B.l) 
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Similarly, the generator L2 of Lorentz boosts in the zj/-plane and the generator 
Lq of Euclidean rotations in the xy-plane, are expressed by 

L2 = s{dt- ^d^ - i>d.4,) + i (1 + e"'* - + 2^^) ds , 
Lo = s {-dt + 4>d.^ + ^Jd^) + 5(1- + - 2^^) 9, . 

Being the generators of the Lie algebra S0i,2 of the Lorentz group, the operators 
Lo, Li, L2 satisfy the commutation relations: 

[Lq, Li] = —L2 , [io, L2] = Li , [ii, L2] = Lo • 

In particular, the generator Lq + L2 = [Li, Lo + L2] is the generator of transla- 
tions of the coordinate s . 

So far, we have been concerned with the case of a single site. To pass to a 
lattice A with many sites, we take the sum 

Ljj) (a = 0,1, 2) 

of differential operators over all sites. 

By construction, the S0i^2 operators L^ — '^La{j) are symmetries of the 
Berezin measure D/^a • Therefore, they give rise to Ward identities: 

= J DfXA La (e-^-^-F) = {LaF - FLaAp,,) [a = 0, 1, 2) , (B.2) 

which hold for any observable F as long as these expectations exist. By com- 
puting the symmetry-breaking terms from the formula LaA^.^ — e La{j) Zj 
one obtains these Ward identities in the more explicit form 

(LiF) = eJIj ((sinhtj - ^e'^s^ _ g'^' Vij?/'j)F) , 
(L2F) =e^^.(c*^s,F) , (LoF)-G. 

The sum rule (|4.6p now follows from the identity for {L2F) by taking F = e*'Si 
and performing the Gaussian integrals over the fields ip, ip, and s. 

Another important consequence results from making the choice F = e*^ sj . 
Since LqF = — sinhtj -t- e*^ (^s^ -I- "t/jji/jj), it follows from {LqF) = that 

(e*j) = (coshtj -I- sinhtj) = (coshi^- -I- e*J (isj -I- ^jtpj}) = {zj) ^ 1 . (B.3) 
The last step, (zj) = 1, is by Proposition [2] of Appendix ICl 



C SUSY Ward identities 



The action function of our H^'^ model has a global symmetry w.r.t. the Lie su- 
peralgebra g :— osp2|2 (fo'" £ ^ 0). As a result, there exist supersymmetric 
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Ward identities for suitable (o5p2|2 invariant) observables. Although such iden- 
tities are standard material from the theory of localization of supersymmetric 
integrals [T7], we nonetheless give their derivation for completeness here, as the 
said identities play a central role in our analysis. The essence of the argument 
can already be understood at the very special example of a lattice A consisting 
of just a single site. For pedagogical reasons we first handle this simple situation 
and then, in a second step, give the generalization to arbitrary lattices. 

The treatment will be most transparent if we do all calculations using the 
coordinates x, y, ^, rj described at the beginning of Section[21 As stated there, for 
our purposes we may view osp2|2 as the space of first-order differential operators 
D with coefficients that are linear in the variables x,y,S,,'i] and the property 

Di/ = 

of annihilating the quadratic polynomial 

H ^ x'^ + + 2^r] . 

Let Q be the distinguished first-order differential operator defined by 

Q ^ xd,^ ~ yd^ + £,da; + rjdy ■ (C.l) 

Clearly Q is odd, converting even coordinate generators a;, y into odd generators 
^, r] and vice versa. Q is also seen to annihilate H, and thus represents an element 
of osp2|2 • Notice that Q squares to 

= xdy - ydx + S,d,j - ijd^ , 

which is a generator from the Lie algebra part 02 ® sp2 of osp2|2 ■ 

Now recall from Section [2] that our Berezin superintegration form is 

D^i= {2tt)-^ dxdy d^d,,^ o (l + H)-^/'^ . 

Lemma 15 The Berezin superintegration form D^i is Q-invariant, i.e., 

[ DnQf = 
for any bounded smooth superfunction f = f{x,y,S,,ri). 

Proof. Since Q is a first-order differential operator, we have from QH = 
that Q{l + H)-^^^ ^ 0. Therefore, D^iQf = Dfi{l + Hy^^Q{l + H)-^/^f and 

/ D^iQf = i2n)-^ f dxdyd^dr, {xdr, - yd^ -f $9, + r,dy) (1 + Hy^'^f . 

The desired result now follows because 9| = 9^ = and the integral over of 
the total derivatives dx{l + H'^)~^^'^f and dy(\ + vanishes. ■ 

An important property of the differential operator Q is that the joint zero 
locus of its coefficients is the origin x ^ y = Q and ^ = 77 = 0. Denoting the 
origin by o we write f{x = 0, y = 0, ^ = 0, 77 = 0) = /(o). 
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Lemma 16 Let f = fix,y,^,r]) be a smooth superfunction which satisfies the 
invariance condition Qf = and decreases sufficiently fast at infinity in order 
for the integral J^2 DjjL f to exist. Then 

[ Dfif = f{o). 



Proof. The idea is to 'deform' the integrand / (without changing the integral) 
by a factor that locaUzes the integral at o. We will do this deformation by 
multiplication with e~'^^ for some positive real parameter t. Thus we are going 
to show that 

''D^lf = J D^ie-^"f, (C.2) 

independent of r > . The desired result will then follow by taking r — > +oo . 

We begin by observing that the localizing function H is Q-exact: it can be 
written as H = QX with 

A ■.= xT]-y^ 

an odd superfunction. Next, using the relation QH — Q^X — we do the 
following calculation: 



e 



- 1 + (e-^t'?^) - 1) = 1 + Q ( A 



g-r(QA) _ I 



QX 



Here the term in parentheses stands for 

g-r(QA) _ I °° 



Inserting this decomposition into the integral we obtain 



Dfie-^"f^ / Dfif+ / DfifQlX 



QX 



Since our integrand / is Q-invariant by assumption {Qf ~ 0), the second integral 
can also be written as 

j D,fQ[^,^^y J D.Q^fX^^yO, 

which vanishes by Lemma [TS] This already proves (|C.2[) . 

To complete the proof, we consider the effect of a scale transformation (j)*: 
X 1-^ xjyjr, y 1-^ y/V^, f/v^j V ^ v/V^- Note that <j)*H — H/t and the 
Berezin superintegration form dxdy d^dri = Dfi o (1 + i?)^/^ is invariant by 0* . 
The statement of the Lemma now results from taking the limit 

/ D^i f = lim f Dtie~^"f= lim / r,{Dfie-^"f) 
= lim [Dii{l + Hf'^{\ + H/T)-^'^e-"4>*f=f{o), 

T^COJ 
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where the last step is done by verifying the normaUzation integral 

and observing that limT-_>oo 0^/ constant function of value /(o). ■ 

We finally turn to the setting of an arbitrary lattice A. We have a first-order 
differential operator Qj for every site j € A and we now take the symmetry 
generator Q to be the sum of all of these: 

By the same argument as before, one sees that -D/^a is Qj-invariant for all j 
and hence Q-invariant. There still exists H = J^jeAi^] +2/1 + '^^jVj) cind 
A = J2jeAi-'^j Vj ~ Uj Cj) with QX = H . Hence we can still localize the integral 
/ -D/XA F for any Q-invariant function F by deforming with e~'^^ and sending 
T oo. Thus we arrive at the following result which, though valid for any choice 
of coordinate system, will be stated in terms of the horospherical coordinates 
tj, Sj, tpj, V'i used in the body of the paper. 

Proposition 2 For any Q-invariant, smooth and integrable function F of the 
lattice variables tj , Sj , -tpj , ipj the integral of F localizes at the zero-field configu- 
ration tj = Sj = ipj = ipj — (for all j G A): 

[ Dti^F^ F{o) . 

J(R2)|A| 

In particular, for the partition function (14. 5p we have 

Z{(3,e) = [ i^/iAe-^-^- = 1 . 

J(R2)|A| 
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